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PLOTSPEC 

A  FORTRAN  77  Program  For  Plotting  Spectral  Data  From 
The  Varian  CARY  2390  UV-VIS-NIR  Spectrophotometer 

ENTRODUCTION 

Spectral  data  acquisition  from  the  CARY  2300  and  2400  series  spectrophotometers  has 
been  implemented  successfully  using  a  Hewlett-Packard  1000  minicomputer  system  using  the 
FORTRAN  77  program  CARYSPEC,  which  is  described  in  a  separate  report.  The  ability  to 
acquire  high  quality  UV-VIS-NIR  spectra  must  be  complemented  with  flexible  data  analysis 
and  plotting  software  to  make  full  use  of  the  spectroscopic  data.  This  report  describes  a  fully 
tested  FORTRAN  77  program  PLOTSPEC  which  fulfills  the  plotting  requirements  of  most 
users  in  a  simple  to  use  menu  driven  environment.  The  program  supports  Absorbance, 
Extinction  Coefficient,  %  Transmission  and  %  Reflectance  ordinates  vs  linear  Wavelength  (nm) 
or  Energy  (cm“*)  abscissa.  Automatic  data  rescaling  is  provided  between  Absorbance  units  and 
%  Transmission  so  that  either  type  of  data  file  from  CARYSPEC  may  be  plotted  in  the 
ordinate  of  choice.  PLOTSPEC  supports  two  concurrent  spectra,  of  up  to  10001  dau  points 
each,  allowing  difference  spectra  to  be  obtained  with  automatic  scaling  for  differences  in 
concentration  and  pathlength.  This  mode  is  only  valid  for  Absorbance  of  Extinction  Coefficient 
ordinate  choices  and  %  Transmission  data  -files  will  be  automatically  rescaled  before  plotting. 
Data  treatment  in  the  %  Reflectance  mode  is  identical  to  %  Transmission  with  only  the 
plotting  label  altered  to  match  the  use  of  this  measurement  technique.  The  abscissa  axis 
scaling  is  bidirectional  to  suit  individual  preferences  for  the  direction  of  increasing  wavelength 
or  energy  units. 

PLOTSPEC  has  been  developed  for  use  with  an  HP  1000  minicomputer  system  running 
the  RTE-6/VM  operating  system  and  Cl  shell.  In  this  environment  the  program  communicates 
with  Hewlett-Packard  model  9872A  and  7550A  digital  plotters  via  the  IEEE-488  and  RS-232C 
I/O  subsytems,  respectively.  These  communication  functions  are  transparent  to  high  level 
languages  on  the  HP  1000  enabling  simple  READ/ WRITE  statements  from  FORTRAN  to 
control  the  external  plotting  operations.  The  program  resides  in  a  single  32K  word  memory 
segment  and  utilizes  79  pages  of  Extended  Memory  Addressing  (EMA)  memory  for  the  large 
data  arrays.  Since  PLOTSPEC  makes  use  of  very  few  special  features  of  the  HP  1000 
computer  system  the  program  could  be  modified  easily  to  run  on  other  host  systems 
supporting  the  FORTRAN  77  language  and  either  the  IEEE-488  or  RS-232C  interface 
standards. 

.Manuscript  approved  August  9,  1988. 
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IMPLEMENTATION 


1.0  Hardware  Interface: 

The  IEEE-488  interface  for  the  HP  1000  system  is  implemented  with  an  HP  59310B 
interface  card  which  utilizes  4  Logical  Unit  (LU)  addresses  in  the  system.  The  LU  addresses 
used  within  the  program  are  system  dependent  and  are  defined  during  the  system  generation. 
The  interface  card  accessed  by  PLOTSPEC  occupies  LU  addresses  31-34,  corresponding  to 
card  addresses  0-3.  Address  0  enables  a  special  addressing  mode  for  sending  low  level  bus 
command  sequences.  Addresses  1-3  are  predefined  automatic  READ/ WRITE  channels  which 
select  device  numbers  1-3  on  the  bus.  PLOTSPEC  uses  LU  33  to  automatically  address  device 
#2  on  the  bus  (HP  9872A  plotter).  To  simplify  alterations  to  the  plotter  address  the  program 
uses  the  INTEGER  variable  PU  for  the  unit  number  in  all  READ/WRITE  operations  with  the 
plotter. 

The  RS-232C  interface  to  the  HP  7550A  plotter  is  implemented  via  a  standard  HP  1000 
8  channel  multiplexer  card  (MUX)  operating  at  9600  baud  with  8  data  bits  and  XON/XOFF 
software  handshaking  protocol.  PLOTSPEC  does  not  send  device  control  commands  to  set 
these  communications  parameters  within  the  HP  7550A  since  these  are  easy  to  select  from  the 
front  panel  control  of  the  plotter.  In  this  application  the  HP  7550A  plotter  is  operated  in  the. 
STANDALONE  mode  with  DIRECT  connection  to  the  MUX  output  of  the  HP  1000. 

1.1  System  Handshaking: 

The  IEEE-488  subsystem  of  the  HP  1000  is  operated  by  the  RTE  driver  DVA37 
configured  for  ASCH  Data  Record  mode,  enabling  automatic  ASCII  NUMERIC  conversion 
on  I/O.  This  mode  sends  and  expects  to  receive  an  End  Of  Record  (EOR)  with  data 
transmission  in  the  form  of  a  Carriage  Retum/Line  Feed  (CR/LF)  sequence  which  matches  the 
requirements  of  Hewlett-Packard  digital  plotters. 

The  RS-232C  subsystem  is  operated  by  the  RTE  driver  DDVOO  which  also  performs 
automatic  ASCII  ♦  NUMERIC  conversions  for  FORTRAN  LO  statements.  This  allows 
PLOTSPEC  to  control  either  plotter  without  regard  for  differences  in  the  handshaking 
protocols.  Character  labelling  mode  on  these  plotter?  requires  a  special  terminat'^r 
and  PLOTSPEC  issues  the  default  Etx  character.  ASCII  3. 
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1.2  HP-GL  Software  Control  Of  HP  Plotters: 


PLOTSPEC  uses  a  subset  of  the  Hewlett-Packard  Graphics  Language  (HP-GL)  to  control 
the  operations  of  the  HP  9872A  and  HP  7550A  digital  plotters.  The  early  generation  of  HP 
plotters  often  only  allow  scaling  and  plotting  commands  to  use  integer  units  which  restricts  the 
utility  of  the  automatic  scaling  command  'SC*  in  these  plotters.  This  lack  of  utility  is  so 
severe  that  the  operating  manual  for  the  HP  9872A  plotter  does  not  even  document  the 
existence  of  the  'SC*  command  in  the  HP-GL  syntax.  However,  PLOTSPEC  utilizes  a  general 
scaling  procedure  with  the  'SC  command  which  functions  on  all  models  of  HP  plotters, 
maintaining  a  0-10000  unit  scale  on  each  axis  at  all  times.  The  'SC*  command  has  the 
advantage  of  automatically  mapping  these  user  units  onto  the  scaling  points  PI  (lower  left) 
and  P2  (upper  right).  In  order  to  map  the  real  user  units  onto  this  10000  x  10000  unit  grid 

the  program  only  needs  to  use  a  multiplier  (XMULT.YMULT)  for  the  data  on  each  axis.  For 

example,  a  %  Transmission  range  of  0-100%  would  require  YMULT  =  100  for  correct  scaling 
of  the  user's  data.  This  technique  maintains  software  resolution  at  1  in  10000  units  regardless 
of  the  range  of  the  user's  MIN, MAX  values. 

PLOTSPEC  also  supports  bidirectional  plotting  on  the  X-axis  even  though  the  'SC 

scaling  command  of  early  HP  plotters  does  not  permit  setting  Xmin  >  Xmax,  i.e.  numbers 
increasing  towards  the  left.  In  order  to  accomplish  this  task  in  a  general  manner  the  program 
uses  two  Xmin, Xmax  ranges,  0  -  +10000  (increasing  to  the  right)  and  -10000  -  0  (increasing 
to  the  left).  The  change  in  sign  has  the  effect  of  reversing  the  direction  of  plotting.  This 
scheme  is  known  to  work  on  HP  9872A,  7225A,  7470A,  7475A  and  7550A  model  plotters. 

The  HP-GL  commands  used  within  PLOTSPEC  are  summarized  below  together  with  their 
syntax  and  parameter  types.  The  type  labels  INTEGER  and  DECIMAL  should  be  interpreted 
as  their  ASCII  representations.  All  data  sent  to  or  received  from  the  plotters  is  in  ASCII 

code.  Therefore,  where  the  command  syntax  below  indicates  a  program  variable,  a  parameter 
value  enclosed  within  literal  string  delimiters  is  equally  applicable.  The  HP  1000  computer 
system  performs  the  conversion  between  numeric  and  ASCII  representations  automatically  on 
I/O  operations  allowing  liberal  mixtures  of  literal  string  and  numeric  variables  to  appear  in 
HP-GL  instructions,  with  a  few  exceptions.  For  example,  the  character  plot  instruction  'CP* 
fails  on  receipt  of  the  value  0.0  from  a  program  variable  but  happily  accepts  the  literal  '0.0'. 
Where  HP-GL  commands  are  known  fail  to  function  correctly  for  valid  ranges  of  parameters 
PLOTSPEC  first  converts  the  to  their  string  representation  before  sending  them  to  the 

plotters. 


1.3  HP-GL  Commands  Summary; 


’AP’  Automatic  Pen  Pickup 

Syntax  =  'AP;' 

’CP’  Character  Plot.. .moves  the  pen  in  units  of  character  width,  height 

Syntax  =  'CP'.X.'.'Y,';' 

DECIMAL  X.Y 

’DF’  Default  Settings... sets  the  plotter  to  default  parameters 

Syntax  =  'DF;' 

’DI’  Direction... sets  the  labeling  direction 

Syntax  =  'DI'.X.'.'.Y,';'  where  X.Y  =  cos(0),sm(0) 

DECIMAL  X.Y  DEFAULT  =  1,0 

’IP’  Input  P1,P2  Scaling  Points. ..sets  the  size  of  the  plot 

Syntax  =  'IP'.Xl  .'.',Y1  ,',',X2.'.',Y2.';' 

where  (Xl,Yl)  =  lower  left  scaling  coordinate  (plotter  units) 
wnere  (X2,Y2)  =  upper  right  scaling  coordinate  (plotter  units) 

INTEGER  X1,Y1,X2,Y2 

’IW’  Input  Window. ..sets  the  size  of  the  plotting  window 

Syntax  =  'IW.aI  .'.'.Yl  ,',',X2,',',Y2,';' 

where  (X1,Y1)  =  lower  left  soft  clip  coordinate  (plotter  units) 
where  (X2,Y2)  =  upper  right  soft  clip  coordinate  (plotter  units) 

INTEGER  X1,Y1,X2,Y2 

’LB’  Label... draws  the  following  literal  characters 

Syntax  =  'LB', string, Etx  where  Etx  =  CHAR(3)  terminates  LB 

ASCII  string 

’LT’  Line  Type.. .solid  or  broken  lines  styles 

Syntax  =  'LT',PN,',',PL,';'  where  PN  =  Pattern  number  (1-6) 

INTEGER  PN,  DECIMAL  PL  where  PL  =  PuvLcin  length  (0-100%) 
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Output  Actual  Pen  Position... sends  the  current  coordinate  position 
Syntax  =  'OA;'  followed  by  READ  (...)  X.Y.Z 
where  (X,Y)  =  current  position  (plotter  units) 
where  Z  =  0  or  1  for  pen  status  (up/down) 

INTEGER  X.Y.Z 

Output  Scaling  Points  PI, P2... sends  the  coordinates  of  plot  limits 

Syntax  =  'OP;'  followed  by  READ  (...)  X1.Y1.X2.Y2 

where  (XI. Yl)  =  lower  left  scaling  point 

where  (X2.Y2)  =  upper  right  scaling  point 

INTEGER  X1.Y1,X2.Y2 

Plot  Absolute... moves  the  pen  to  an  absolute  coordinate 

Syntax  =  'PA'.X.'.'.Y,';'  where  X.Y  =  coordinates  (plotter  units) 

HP  9872A;  INTEGER  X.Y 

HP  7550A;  INTEGER  X.Y  or  DECIMAL  X.Y  if  user  scaling  is  on 

Pen  Down... sets  the  pen  on  the  plotting  surface 

Syntax  =  'PD;' 

Pen  Up. ..picks  up  the  pen  from  the  plotting  surface 

Synux  =  'PU;' 

Scale  Plot.. .maps  user  units  for  X.Y  axes  onto  P1,P2  scaling  points 

Syntax  =  'SC'.X1.'.'.X2.'.',Y1.'.'.Y2.';' 

where  (Xl.Yl)  =  coordinate  of  scaling  point  PI  (user  units) 

where  (X2.Y2)  =  coordinate  of  scaling  point  P2  (user  units) 

HP  9872A:  INTEGER  X1.Y1.X2.Y2 
HP  7550A;  DECIMAL  X1.Y1.X2.Y2 

Size  Relative... sets  the  character  size  in  proportion  to  axes  scaling 

Synux  =  'SR', WIDTH.'. '.HEIGHT.';' 

DECIMAL  WIDTH.HEIGHT 


’SP’  Select  Pen.. .picks  up  the  pen  from  selected  stall  or  carousel  position 

Syntax  =  'SP*,X,';'  where  X  =  pen  number  (0-8) 

LNTEGER  X 

’TL’  Tick  Length... sets  the  length  of  positive  tick  marks  for  axes 

Syntax  =  'TL',X,';'  where  X  =  %  of  full  scale  (0-100) 

DECIMAL  X 

’VS’  Velocity  Select.. ..sets  the  pen  down  speed 

Syntax  =  'VS’.X, where  X  =  velocity  (cm/sec) 

INTEGER  X 

’XT’  X- tick... draws  a  vertical  tick  mark  at  the  current  pen  position 

Syntax  =  'XT;' 

’YT’  Y-tick.. .draws  a  horizontal  tick  mark  at  the  current  pen  position 

Syntax  =  'YT 

1.4  Default  Plotting  Parameters: 

PLOTSPEC  initializes  the  HP  plotters  using  the  'DF'  default  parameters  command  and 

then  sets  up  the  following  plotting  style: 

LABELS  The  initial  label  direction  is  set  to  horizontal  with  nominal  character  fields  of 
1x2  cm  and  the  size  relative  option  is  selected. 

PEN  Pen  #1  is  selected  with  automatic  pen  pickup  and  set  to  the  LT  state.  The 

pen  down  velocity  is  set  to  5  cm/sec  for  best  line  quality  with  felt  tip  pens 
on  graphics  paper. 

SCALE  Initial  scaling  coordinates  for  the  HP  y872A  plotter  are  (1350,1000)  and 
(8500,7000)  corresponding  to  the  NOTEBOOK  size  plot  on  A-size  paper 
while  the  plotting  window  is  set  to  the  hard  clip  limits. 


TICKS 


Tick  marks  are  set  to  +1%  of  full  scale. 


2.0  Purpose  Of  PLOTSPEC: 


This  program  is  intended  to  provide  a  flexible  plotting  facility  for  spectral  data  acquired 
from  the  Cary  2390  spectrophotometer  using  the  program  CARYSPEC,  which  is  described  in  a  ■ 

separate  report.  The  HP  1000  computer  system  provides  a  large  EMA  area  to  hold  data 
arrays  and  PLOTSPEC  makes  use  of  this  feature  to  access  two  concurrent  spectra  of  up  to 
10001  points  each.  This  allow  for  difference  spectra  to  be  plotted  with  relative  ease  compared 
with  the  experimental  difficulties  incurred  with  real  time  subtraction.  While  ihe  raw  data  from  I 

the  spectrophotometer  consists  of  Absorbance  or  %  Transmission  measurements  vs  Wavelength 
(nm),  PLOTSPEC  allows  for  presentation  of  the  data  in  more  meaningful  units  as  well.  The 
abscissa  may  be  scaled  in  Wavelength  (nra)  or  Wavenumbers  x  10"^  with  increasing  values 
towards  the  left  or  right.  The  ordinate  may  be  presented  in  Absorbance,  Extinction  Coefficient 
(Molar  Absorptivity),  %  Transmission  or  %  Reflectance.  PLOTSPEC  automatically  rescales  the 
original  ordinate  data  to  match  the  choice  of  plotting  ordinate.  Furthermore,  automatic  range, 
tick  and  label  position  selection  is  provided  for  fast  and  efficient  setup  of  a  neat  and  tidy 
plot.  However,  the  user  can  easily  alter  the  default  selections  in  the  menu  driven  plotting 
environment. 

PLOTSPEC  also  provides  some  useful  annotation  facilities  for  documenting  the  plot  and 
for  producing  presentation  material.  A  digitizing  routine  is  interwoven  with  the  annotation 
commands  to  enable  semi-automatic  labelling  of  peak  maxima  with  the  ease  of  a  point  and 
shoot  technique.  More  accurate  determinations  of  peak  maxima  can  also  be  made  using  a 
page  scrolling  data  listing  routine,  though  the  digitizing  facility  seems  to  be  perfectly  adequate 
in  most  cases.  The  default  plotting  mode  produces  a  small  .NOTEBOOK  size  plot  on  A-size 
paper,  together  with  a  list  of  spectral  parameters  for  documentation  purposes.  A  FULL  size 
plotting  mode  can  also  be  selected  to  produce  a  larger  plot  on  A-size  paper  without  a 
parameter  list.  Both  modes  include  sufficient  space  for  a  centered  title  over  the  plot. 

Additionally,  the  program  provides  a  USER  SET  plotting  size  with  default  boundaries  suited 
for  a  full  size  plot  on  B-size  paper.  This  mode  is  ideal  for  producing  poster  presentation 
material.  PLOTSPEC  supports  output  to  two  generic  types  of  Hewlett-Packard  digital  plotters, 
the  older  HP  9872A  and  the  more  recent  HP  7550A  model.  However,  PLOTSPEC  is  designed 
to  function  with  all  generations  of  HP  plotters  with  only  minor  differences  in  the  position  of 
the  plotting  origin. 
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2.1  Program  Structure: 


PLOTSPEC  comprises  a  large  main  program  unit  containing  the  console  menu  displays, 
string  data  for  the  instrument  settings  variables  and  a  number  of  subroutines  that  perform 
string  processing,  input  validation  and  communication  with  the  HP  9872A  and  HP  7550A 
digital  plotters.  The  main  program  unit  is  responsible  for  all  the  logic  flow  and  the 
subroutines  execute  specific  support  tasks,  which  are  summarized  below; 

The  main  program  unit  of  PLOTSPEC  comprises  6  distinct  segments  of  code  to  carry 
out  the  the  functions  of  disk  file  data  retrieval,  instrument  settings  display,  plotting  mode 
setup  and  actual  plotting  operations.  The  code  fragments  appear  under  the  following  assigned 
labels:  MENU,  READ,  SETTINGS,  DIFF,  PLOT  and  EXIT. 

MAIN  PROGRAM 


MENU; 

This  is  the  first  and  main  control  menu  of  the  program,  selecting  entry  to  the  data 

retrieval,  instrument  settings,  plotting  mode  and  exit  routines.  The  choices  are  as  follows: 

’Rl’  .....Read  Spectrum  #1 

This  command  branches  to  label  READ  and  the  program  performs  logical  tests  for 
the  presence  of  a  valid  spectrum  in  memory  before  allowing  previous  data  to  be 
overwritten  by  a  new  file.  The  data  file  format  is  listed  below  in  Table  I.  Files 

are  read  into  the  first  dimension  of  the  data  arrays  and  may  be  plotted  using  the 
'PI '  command. 

’R2’  . Read  Spectrum  #2 

This  command  performs  the  same  operation  as  'Rl '  but  the  file  is  read  into  the 
second  dimension  of  the  data  arrays.  The  second  spectrum  can  be  plotted 

independently  of  the  first  or  used  for  difference  spectroscopy. 

’IS’  . Instrument  Settings 

This  option  branches  to  label  SETTINGS  and  prompts  for  the  number  of  the 
spectrum.  This  routine  allows  for  display  of  the  most  important  operating  conditions 
of  the  Cary  2390  spectrophotometer  during  acquisition  of  the  selected  spectrum. 
This  provides  an  on-line  reference  when  comparing  a  number  of  spectra. 


’PS’  . Plotter  Selection 

This  option  allows  the  user  to  select  either  the  HP  9872A  (default)  or  HP  7550A 
plotter  as  the  outpui  device.  Digitizing  is  more  convenient  with  the  older  HP 
9872A  plotter  but  output  quality  and  plotting  speed  are  much  improved  with  the 
HP  7550A. 

’PI’  . Plot  Spectrum  #1 

Selection  of  this  command  branches  to  label  PLOT  which  performs  a  check  for  the 
presence  of  a  valid  file  in  memory.  The  routine  then  starts  the  plotting  mode 
setup  for  the  first  spectrum  with  selection  of  ‘he  abscissa  and  ordinate  scaling 
modes.  The  user  has  a  choice  of  Nanometres  and  Wavenumbers  for  the  abscissa 
and  a  choice  of  ABSORBANCE,  Extinction  Coefficient.  %  Transmission  and  % 
Reflectance  for  the  ordinate.  The  program  then  proceeds  to  scale  the  data  and  set 
default  plotter  operating  modes  before  presenting  a  menu  of  plotting  options.  The 
original  data  file  noay  have  been  recorded  in  either  Absorbance  or  %  Transmission 
units  so  the  data  will  be  rescaled  if  the  choice  of  plotting  ordinate  differs  from  the 
file  ordinate. 

’P2’  . Plot  Spectrum  #2 

This  command  performs  the  same  operations  as  'PI  ’  but  uses  data  residing  in  ih . 
second  dimension  of  the  data  arrays. 

’DS’  . Difference  Spectrum  (1-2) 

Selection  of  this  option  causes  a  branch  to  label  DIFF  which  first  checks  for  the 
presence  of  two  valid  spectra  in  memory  and  then  tests  for  matching  of  their 
wavelength  range  and  step  size  (nm).  If  the  spectra  are  matched  the  prog’-am 
continues  to  label  PLOT  in  the  normal  fashion.  In  this  case,  however,  all  plotting 
operations  will  use  the  scaled  difference  between  the  two  data  sets.  The  data  from 
the  second  spectrum  are  scaled  for  differences  in  concentration  or  pathlength 
before  subtraction  during  plotting. 

’EX’  . Exit 

The  final  option  causes  an  unconditional  branch  to  the  label  EXIT  which 
terminates  the  program. 
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TABLE  I 


Data  File  Format 


Line  File  Variables 


Format  Type 


1 

2 

3 

4 


5 

6-54 

55 


56-/ 


/-eof 


TITLE(2) 

DATE(2) 

FMIN(2),FMAX(2),FSTEP(2).CONC(2), 

PATH(2) 

ORD(2).ABSC(2).CELL(2).CYCLE(2), 

SAMPLE(2).WAVE(2),TIMER(2). 

TEMP(2),DIST(2) 

J.K,NDATA(2,I) 

PARAM(2,I) 

VARIABLE(2.I) 

A(2.I) 

W(2.I) 


CHARACTER  (A72) 
CHARACTER  /A8) 
REAL  (•) 

REAL  (•) 


INTEGER  (13,13,16) 
INTEGER  (12) 
REAL  (•) 

REAL  (•) 

REAL  (•) 


I 


a:  Disk  data  files  are  read  into  either  half  of  the  data  arrays 


b:  (•)  indicates  free  field  format 


READ; 


This  section  of  PLOTSPEC  opens  a  disk  file  previously  stored  by  the  data  acquisition 
program  CARYSPEC.  The  data  are  stored  in  ASCII  code  and  contain  a  complete 
description  of  the  instrument  parameter  settings  as  well  as  the  spectral  data.  The  format 
of  the  data  file  is  listed  above  in  Table  I.  The  program  prompts  the  user  for  both  the 
filename  and  subdirectory.  The  latter  defaults  to  the  user's  directory  if  only  a  RETURN 
character  is  entered.  If  a  directory  name  is  given  the  program  builds  a  complete  Cl 
pathname  for  use  in  the  OPEN  statement.  If  the  file  cannot  be  found  or  the  file  is 
already  open  the  program  displays  an  error  message  and  then  returns  to  the  main 
MENU.  Otherwise,  the  file  is  read  and  the  file  descriptors  are  displayed  on  the  console 
while  the  remainder  of  the  data  are  being  transferred.  After  completion  of  the  data 
transfer  the  program  prompts  for  changes  to  the  concentration  or  pathlength  variables. 
This  feature  allows  for  on-line  corrections  to  the  file  variables  before  plotting,  though 
permanent  corrections  should  be  made  using  the  file  editing  program  EDITSPEC.  The 
program  then  returns  to  the  main  MENU. 


SETTINGS: 

This  section  of  code  displays  a  list  of  the  most  important  instrument  settings  of  the  Cary 
2390  spectrophotometer  during  acquisition  of  the  specified  spectrum.  The  program 
includes  a  large  amount  of  string  data  for  the  various  settings  in  the  CHARACTER  array 
Pstr.  The  INTEGER  array  PARAM  is  used  as  an  index  to  this  string  data  while  numeric 
data  for  the  table  are  obtained  from  the  REAL  array  VARIABLES. 


DIFF; 

I  Entry  into  this  section  of  code  is  made  prior  to  plotting  in  order  to  verify  that  the  data 

in  memory  are  suitable  for  difference  spectral  plots.  The  program  first  checks  that  there 
are  two  spectra  present  and  proceeds  to  test  for  exact  matching  of  the  wavelength  limits 
and  step  size.  The  editing  program  EDITSPEC  includes  routines  for  preparing  subset  data 
I  files  to  support  the  difference  plotting  mode  of  PLOTSPEC. 
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PLOT: 


This  section  of  PLOTSPEC  allows  the  user  to  set  up  the  axes  scaling  to  suit  the  type  of 
data  to  be  plotted.  If  the  selected  spectrum  is  not  present  in  memory  an  error  message 
is  displayed  and  the  program  returns  to  the  main  MENU.  Otherwise,  a  prompt  is 
displayed  for  the  choice  of  abscissa.  PLOTSPEC  supports  two  abscissa  modes.  Wavelength 
(nm)  and  Wavenumbers  x  10~3.  After  this  selection  the  program  sets  a  number  of  X-axis 
plotting  variables  and  calls  Subroutine  Xaxis  to  setup  the  default  tick  and  label 
positioning.  The  abscissa  is  bidirectional  allowing  the  user  to  select  the  direction  of 
increasing  values.  However,  the  default  direction  is  for  values  to  increase  towards  the 
right. 

The  following  prompt  is  for  the  choice  of  ordinate  mode.  Data  may  be  plotted  in  a 
choice  of  ABSORBANCE,  Extinction  Coefficient,  %  Transmission  and  %  Reflectance. 
The  latter  pair  are  identical  e  rept  for  the  label  on  the  Y-axis.  Data  files  from  the 
CARYSPEC  program  may  contain  either  Absorbance  or  %  Transmission  values  so  that 
rescaling  is  provided  automatically  to  suit  the  user’s  choice  of  plotting  ordinate.  If  the 
Difference  Spectrum  mode  has  been  invoked  only  the  Absorbance  and  Extinction 
Coefficient  ordinates  are  valid.  If  the  plotting  ordinate  uses  Absorbance  data  the  program 
will  search  the  data  array  for  the  maximum  value  to  provide  automatic  ranging  on  the 
Y-axis.  If  Extinction  Coefficient  mode  is  specified  the  program  also  calculates  a  suitable 
Y-axis  exponent.  Otherwise,  the  Y-axis  defaults  to  0-100%  limits.  The  default  tick  and 
label  positions  are  then  set  by  calling  Subroutine  Yaxis. 

The  user  l»  then  prompted  to  check  that  the  selected  plotter  is  connected  and  turned  on 
before  proceeding  to  send  plotter  commands.  Finally,  the  user  is  presented  with  the 
plotting  menu  with  default  settings  for  axis  and  label  sizes,  line  type,  pen  velocity,  the 
location  of  tick  marks  and  axis  labels  and  the  command  interface  for  plotting,  annotating 
and  digitizing. 


EXIT: 

The  final  ponion  of  PLOTSPEC  issue  an  erase  command  line  instruction  and  then 
terminates  the  program  unconditionally. 


2.2  Plotting  Menu  Commands: 


The  interactive  plotting  control  menu  is  designed  both  for  flexibility  and  ease  of  use  in 
setting  quite  a  large  number  of  plotting  parameters.  On  most  occasions,  the  default  limits  and 
spacing  between  tick  marks  and  labels  will  be  perfectly  satisfactory.  However,  tne  user  has 
complete  control  over  these  parameters  to  suit  individual  preferences.  The  various  command 
selections  are  listed  below; 

’PL’  . Plot  Limits 

This  command  sets  the  actual  range  of  the  data  file  to  be  plotted.  The  default 
value  is  the  full  scan  range.  If  the  plotting  scale  is  made  smaller  than  the  scan 
range  the  plot  limits  are  reduced  to  match  in  order  to  prevent  attempts  to  plot  off 
scale  data.  If  the  plotting  scale  is  made  larger  than  the  scan  range  the  plot  limits 
are  set  to  the  actual  scan  limits. 

’XS’  . X-scale 

The  range  of  the  X-axis  plotting  scale  defaults  to  the  scan  range  of  the  spectrum 
but  may  be  made  larger  or  smaller  than  the  scan  limits.  The  XMIN.XMAX  values 
may  be  entered  in  reverse  order  to  produce  a  plot  with  values  incre.asing  towards 
the  left.  These  limits  are  passed  to  Subroutine  Xaxis  to  reset  the  plot  limits,  tick 
marks  and  label  positions  automatically. 

’XT’  . X-tick  Spacing 

This  option  allows  for  customizing  the  X-axis  tick  mark  spacing  with  prompts  for 
First,  Last  and  Space.  The  limits  are  validated  to  lie  within  the  current  X-scale 
range  and  the  label  positions  are  automatically  reset  to  match  the  tick  positions. 

’XL’  . X- label  Spacing 

This  option  allows  the  user  to  reset  the  automatic  label  positioning.  The  entries  for 
First,  Last  and  Space  are  treated  in  a  similar  manner  to  the  tick  marks. 

’YS’  . Y-scale 

The  range  of  the  Y-axis  plotting  scale  defaults  to  the  MIN, MAX  scale  appropriate 
for  the  selected  ordinate.  For  %R  and  %T  scales  the  default  scale  is  0-\00%.  For 
Absorbance  and  Extinction  Coefficient  ordinates  autoranging  is  used  to  fit  the  entire 
spectrum  on  scale.  The  scale  range  is  passed  to  Subroutine  Yaxis  to  reset  the  plot 
limits,  tick  marks  and  label  positions  automatically. 
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’YT’  . Y-tick  Spacing 

This  option  allows  for  customizing  the  Y-axis  tick  mark  spacing  with  prompts  for 
First,  Last  and  Space.  The  label  positions  are  automatically  reset  to  match  the  tick 
marks.  No  input  validation  is  provided  for  the  limits  as  it  is  assumed  the  user  can 
read  the  updated  console  display  before  plotting. 

’YL’  . Y- label  Spacing 

This  option  allows  the  user  to  reset  the  automatic  label  positioning.  The  entries  for 
First,  Last  and  Space  are  treated  in  a  similar  manner  to  the  tick  marks. 

’CS’  . Char.Size 

The  size  of  character  labels  can  be  altered  to  suit  individual  preferences.  However, 
the  size  should  not  be  made  much  larger  than  the  default  fields  (1x2  cm)  when 
using  the  NOTEBOOK  or  FULL  size  axes  definitions  since  the  labels  may  not  fit 
within  the  hard  clip  limits  of  the  plotters. 

TV’  . Pen  Velocity 

Plotting  quality  is  very  dependent  on  the  pen  speed,  the  type  of  pen  and  plotting 
medium  chosen.  The  program  uses  a  default  pen  speed  of  5  cm/sec  for  both  the 
axes  and  spectrum  plots.  The  user  can  alter  these  within  the  range  1-36  cm/ sec. 
While  the  manual  for  the  HP  7550A  plotter  recommends  a  speed  of  50  cm/sec  for 
the  combination  of  felt  tip  pen  and  graphics  paper  there  is  no  doubt  that  this  is 
determined  more  by  the  quest  for  throughput  than  quality.  Roller  ball  pens  also 
work  well  at  the  5  cm/sec  pen  speed.  Thus,  there  is  little  need  to  alter  the 
default  values  in  most  cases  though  transparency  plots  can  benefit  from  a  higher 
speed  (10  cm/sec). 

’LT’  . Line  Type 

The  plotters  have  6  inbuilt  line  type  functions.  The  default  value  specifies  a 
SOLID  line.  Since  PLOTSPEC  can  plot  multiple  spectra  on  the  same  graph  the 
BROKEN  line  type  function  is  quite  useful.  The  various  patterns  can  be  found  by 
reference  to  the  plotter  manuals.  The  length  of  the  patterns  is  entered  as  a  %  of 
full  scale  and  a  useful  guide  is  to  select  the  same  length  as  the  pattern  number. 
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’AX’  . Axes  Size 


’AN’ 


’PA’ 


’PS’ 


’DD 


The  program  includes  three,  predefined  plot  types  NOTEBOOK,  FULL  and  USER 
SET.  The  plotting  boundaries  vary  to  suit  different  functions  and  paper  sizes.  The 
default  NOTEBOOK  size  is  a  small,  A-size  plot  which  allows  room  for  labelling  a 
list  of  spectral  file  parameters  on  the  right  hand  margin.  The  FULL  size  plot 
nicely  fills  the  plotting  area  of  A-size  paper  with  sufficient  space  for  a  centered 
Title  over  the  plot.  The  USER  SET  size  allows  for  custom  size  plots  up  the  the 
limits  of  B-size  paper.  This  mode  has  predefined  limits  which  produce  a  FULL 
size  plot  on  B-size  paper  but  the  scaling  coordinates  can  be  easily  altered  by  the 
user  to  suit  individual  needs. 

. Annotate  Plot 

This  command  allows  the  use  to  enter,  move  and  plot  Labels,  a  list  of  spectral 
Parameters  or  a  centered  Title.  The  Title  defaults  to  the  data  file  descriptor  but 
the  user  can  enter  a  new  Title  if  required.  The  Parameters  option  is  only  valid  for 
the  NOTEBOOK  size  plot  which  allocates  room  for  these  descriptors  in  the  right 
hand  margin.  The  Labelling  options  allows  for  entry  of  descriptive  text  (up  to  40 
chars.)  which  may  be  positioned  and  drawn  anywhere  on  the  plotting  surface.  The 
character  size  and  labelling  direction  are  selectable  while  the  plotting  mode  supports 
centered,  left  and  right  justification.  The  character  size  for  labels  defaults  to 
0.6* Normal  plot  size. 

. Plot  Axes 

This  command  draws  the  axes  grid  according  to  the  selected  tick  and  label 
positions  listed  in  the  plotting  menu. 

. Plot  Spectrum 

Before  proceding  with  the  plot  this  command  pauses  for  the  possible  entry  of  a 
step  multiplier.  This  ailows  the  user  to  speed  plotting  for  very  large  data  files. 

. Digitize/Display  Data 

The  program  supports  manual  use  of  the  plotters  for  digitizing  peak  positions  from 
the  plot.  This  routine  will  allow  the  user  to  annotate  the  plot  at  these  positions 
with  either  the  X  coordinate,  the  Y  coordinate  or  Both.  Labels  are  drawn  with 
small  characters  in  the  vertical  direction  and  may  be  plotted  in  centered,  left  or 
right  justified  format.  In  this  case,  the  annotation  routine  provides  offsets  to  avoid 
writing  over  the  spectrum. 
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2.3  COMMON  DATA: 


All  COMMON  variables  used  by  PLOTSPEC  are  held  in  named  COMMON  blocks.  The 
variables  contained  within  the  COMMON  blocks  are  listed  below; 

/DATA/  Contains  EMA  REAL  Arrays  of  spectroscopic  data 

A  Array  containing  the  Absorbance  or  %T  data 

Values  read  from  disk  data  file 
Values  used  in  main  program  unit 

W  Array  containing  the  wavelength  data  (nm) 

Values  read  from  disk  data  file 
Values  used  in  main  program  unit 


/DIGIT/  Contains  REAL  and  INTEGER  variables  for  scaling  digitized  data 

FACTOR  REAL  variable  for  scaling  absorbance  to  extinction  coefficient  units 
Values  set  in  main  program  unit  (from  dau  file  variables) 

Values  used  in  main  program  unit  and  Subroutine  Digitize 
XMIN.  REAL  variables  containing  the  current  X-axis  limits 
XMAX  Values  set  in  main  program  unit 

Values  tised  in  main  program  unit.  Subroutine  Digitize,  Subroutine  Order  and 
Subroutine  Xaxis 

YOFF,  REAL  variables  containing  the  current  Y-axis  offset  &  scale  length 
YSCALE  Values  set  in  main  program  unit 

Values  used  in  main  program  unit.  Subroutine  Digitize,  Subroutine  Ore''?  and 
Subroutine  Yaxis 

IX,IY  INTEGER  variables  containing  the  coordinates  of  the  digitized  point 

XI, Y1  INTEGER  variables  containing  the  coordinates  of  the  scaling  points  P1,P2 

X2,Y2  Values  set  in  main  program  unit  and  Subroutine  Digitize 

Values  used  in  main  program  unit  and  Subroutine  Digitize 
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/XPARAM/  Contains  REAL  variables  for  plotting  control 


LMIN,  Current  X-axis  plotting  limits  for  validation  of  Tick  or  Label  entries 
LMAX 

WN.WX  Current  plotting  limits  for  Wavelength  or  Wavenumber  scale 

Xtick,  X-axis  tick  mark  spacing 

FXT.LXT  First  and  Last  X-axis  tick  marks 

Xlabel,  X-axis  label  spacing 

FXL.LXL  First  and  Last  X-axis  labels 

Values  set  in  Subroutine  Xaxis 


Values  used  in  main  program  unit 


/YPARAM/  Contains  REAL  variables  for  plotting  control 

Ytick,  Y-axis  tick  mark  spacing 

FYT.LYT  First  and  Last  Y-axis  tick  marks 
Ylabel  Y-axis  label  spacing 

FYL.LYL  First  and  Last  Y-axis  labels 

Values  set  in  Subroutine  Yaxis 
Values  used  in  main  program  unit 
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2.4  SUBROUTINES: 


PLOTSPEC  uses  subroutines  to  perform  specific  tasks  which  are  required  more  than 
once,  including  string  manipulations,  input  validation,  digitizing  and  data  rescaling.  The 
purpose  and  calling  sequences  are  listed  below: 


Center(TITLE) 

Prints  a  string  on  the  user  console  centred  within  a  72  column  line. 

TITLE  CHARACTER*72  string,  contents  set  by  calling  unit 

CALLED  BY:  Main  program  unit 
CALLS :  None 


hCoeff  (E,String,K) 

Converts  a  numeric  extinction  coefficient  into  a  rounded  string  with 
1  exponent,  if  required.  Special  handling  is  executed  for  0.0000  and  1 .0000 

I  values. 

u 

j  E  REAL  variable  input  from  calling  unit 

t 

i 

String  CHARACTER* (•)  output  containing  floating  point  string 

^  DIMENSION  is  set  by  the  calling  unit 

I  K  INTEGER  variable  input  from  calling  unit  specifying  the  number  of  digits  of 

'  precision  required  in  string.  Output  value  returns  the  number  of  characters  in 

i 

string  to  the  calling  unit. 

I 

I  CALLED  3Y:  SUBROUTINE  Digitize 

!  CALLS;  SUBROUTINE  Str 
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Digitize(X,Xcode,Y,Ext,K) 

Reads  the  current  pen  position  from  the  plotter  and  converts  the  coordinates 
to  the  user's  units.  The  routine  also  converts  the  Y-coordinate  to  extinction 
coefficient  format.  The  main  program  determines  whether  the  latter  value  is 
valid  for  the  ordinate  mode  in  use. 

X  REAL  variable  output  containing  the  abscissa  value  in  user  units 

Xcode  CHARACTER  variable  input  used  to  select  scaling  of  the  abscissa  units 

Y  REAL  variable  output  containing  the  ordinate  value  in  user  units 

Ext  CHARACTER* (•)  output  string  containing  the  extinction  coefficient 

DIMENSION  is  set  by  the  calling  unit 

K  INTEGER  variable  input  to  specify  the  number  of  digits  precision  in  Ext 

Output  value  contains  the  number  of  characters  in  the  returned  string 

CALLED  BY;  Main  program  unit 

CALLS:  SUBROUTINE  Coeff 


». 

Exponent(N,Mult) 

Converts  the  scale  factor  (power  of  10)  for  extinction  coefficient  mode  into 
string  exponents  for  labelling  the  Y-axis  in  the  screen  menu  and  on  the 
plotted  graph  {e.g.  'lE-3'  or  '1  ...’). 

N  REAL  variable  input,  power  of  10 

Mult  CHARACTER* (*)  output  string  containing  multiplier 

DIMENSION  is  set  by  the  calling  unit 

CALLED  BY;  Main  program  unit 

CALLS:  SUBROUTINE  Str  • 


I 
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Lme(N) 

Prints  a  line  of  characters  to  the  user  console  N  columns  wide  and 
centred  within  a  72  column  line. 

N  INTEGER  variable  input  from  calling  unit 

CALLED  BY;  Main  program  unit 
CALLS :  None 


Order(MIN,MAX) 

Validates  the  input  arguments  so  that  MIN  <  MAX 

MIN, MAX  REAL  variables  input  and  output  with  order  swapped  if  required 

CALLED  BY:  Main  program  unit 
CALLS:  None 


Rescale(J,K,ND) 

Converts  ordinate  data  array  between  Absorbance  units  and  %  Transmission. 
The  data  array  A(J,I)  is  accessed  from  EMA  COMMON. 

J  INTEGER  variable  input  (1  or  2)  to  specify  the  spectrum 

K  INTEGER  variable  input  (0  or  1)  to  specify  the  scaling  procedure.  The 

original  value  for  K  is  obtained  from  PAIL*.M(J,1)  in  the  data  file.  On 
return  to  the  main  program  unL  the  latter  parameter  is  altered  to  reflect  the 
current  data  format. 

ND  INTEGER  variable  input  containing  the  number  of  data  points  in  AiJ.I) 


CALLED  BY:  Main  program  unit 

CALLS:  None 
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Str(VALUE,String,PREC) 

Performs  a  conversion  frori  numeric  value  to  a  string  number  for  floating 
point  numbers  only  with  up  to  12  digits  precision. 

VALUE  REAL  variable  input  to  be  processed  by  the  routine 

String  CHARACTER*14  string  output  corresponding  to  VALUE 

PREC  INTEGER  variable  input  to  set  the  rounding  precision  for  string 

CALLED  BY:  Main  program  unit,  SUBROUTINE  Coeff,  SUBROUTINE  Exponent 

CALLS:  None 

Upper(Code) 

Performs  a  check  for  lower  case  characters  in  a  string  of  arbitrary  length 
and  converts  to  upper  case  if  necessary. 

Code  CHARACTER*(*)  variable  passed  into  routine  and  UPPER  case  on  exit 

DIMENSION  is  set  by  the  calling  unit 

CALLED  BY:  Main  program  unit 

CALLS:  None 

EXTENSION:  LEN(5iring)  function,  an  HP  extension  to  FORTRAN  77 

Val(String,  VALUE) 

Performs  a  conversion  from  string  to  numeric  value  for  a  string  number 
containing  up  to  10  digits.  This  is  more  than  required  by  PLOTSPEC. 

String  CHARACTER* (•)  string  input  to  be  processed  by  routine 

DIMENSION  is  set  by  calling  unit 

value  REAL  variable  output 

CALLED  BY:  Main  program  unit 

CALLS:  None 
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Wait(DELAY) 


Performs  a  loop  which  tests  the  system  clock  until  DELAY  seconds  have 

elapsed.  The  routine  does  not  make  provision  for  the  special  case  at  the 

I  transition  to  2400  hours. 

DELAY  REAL  variable  holding  the  value  of  the  delay  period  in  seconds 

I  CALLED  BY:  Main  program  unit 

CALLS:  FUNCTION  Time(I) 

I 

Xaxis(XL,XH) 

Performs  both  limits  settings  for  the  X-axis  range  variables  and  selection  of 
the  default  Tick  and  Label  positions.  The  routine  starts  with  Xtick  set  to  i/io 
of  the  full  scale  range  and  then  compares  this  with  a  table  of  preferred 

values.  If  the  range  is  anomalously  small  or  large  the  Tick  and  Label 
positions  are  based  on  the  i/io  scale  spacing.  Otherwise,  a  neat  and  tidy 
spacing  is  selected  from  the  table  values  with  a  label  at  every  other  tick 
position. 


Yaxis(YSCALE) 

Performs  Tick  and  Label  position  selection  using  preferred  table  values,  if 
possible.  The  routine  start  with  Ytick  equal  to  i/io  of  YSCALE.  If  the 
ordinate  uses  Extinction  Coefficient  scaling  or  the  range  is  anomalously  small 
or  large  the  Tick  and  Label  positions  are  based  on  the  i/io  scale  spacing. 
Otherwise,  a  neat  and  tidy  spacing  is  selected  from  the  table  values  'vith  a 
label  at  every  other  tick  position. 


O'* 


2.5  FUNCTIONS: 


PLOTSPEC  uses  only  one  function  subprogram  that  makes  an  EXEC  call  to  read  the 
system  time. 

Time(I) 

Performs  an  EXEC  call  to  read  the  system  clock  and  converts  the  reading  to 
seconds  and  centiseconds. 

I  Dummy  argument 

CALLED  BY;  SUBROUTINE  Wait  only 

CALLS;  EXEC(ICODE,ITIME)  system  level  command 

PROGRAM  CODE 

3.0  Source  Code  Availability: 

The  source  code  for  program  PLOTSPEC  is  a  73K  ASCII  text  file  available  on  either  a 
Hewlett-Packard  cartridge,  9  track  tape  or  an  IBM  360K  format  floppy  disk.  All  requests 
should  be  accompanied  by  the  blank  medium  desired.  A  printed  copy  of  the  source  code  is 
listed  below. 

3.1  Variable  Names  And  Usage: 

A  complete  listing  of  the  INTEGER,  REAL,  REAL  Array  and  CHARACTER  variables 
for  the  MAIN  segment  of  PLOTSPEC  is  given  below  in  Tables  II,  III,  IV  &  V,  respectively. 
The  subroutines  use  the  same  names  as  the  main  program  for  the  same  variables.  Additional 
variables  in  the  subroutines  and  simple  integers,  I-N,  are  not  documented  since  their  usage  is 
rather  obvious.  The  logical  variable  MATCH  is  used  within  the  program  when  comparing  two 
spectra  for  the  difference  spectrum  mode. 
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Table  II 


Glossary  of  INTEGER  Variables 


Name 

Description 

Value 

DIFF 

Assigned  Label  -  Difference  Spectrum 

370 

EXIT  • 

Assigned  Label  -  Terminate  Program 

9000 

MENU 

Assigned  Label  -  Main  Control  Menu 

10 

PLOT 

Assigned  Label  -  Plot  Spectrum,  1  or  2 

400 

READ 

Assigned  Label  -  Read  Spectrum,  1  or  2 

570 

SETTINGS 

Assigned  Label  -  Scan  Conditions,  1  or  2 

300 

Ascn 

ASCn  equivalent  of  digits  in  Str 

48-57 

FINISH 

Index  of  final  abscissa  value  to  plot 

1-10001 

IX.IY.IZ 

Plotter  position  variables  READ/WRITE 

0-16000 

NCOL 

Number  of  screen  columns  in  menu  display 

50-70 

ND 

Number  of  data  points  in  plotted  spectrum 

1-10001 

NP 

Number  of  parameters  to  read  from  file 

49 

NV 

Number  of  variables  to  read  from  file 

14 

PREC 

Precision  for  rounding  function  in  Str 

3.4 

PN 

Plotter  model  selection  number 

1,2 

PU 

Plotter  Logical  Unit  #,  READ/WRTTE 

17,33 

PXOFF 

Plotter  X-zero  position  offset 

0,200 

PYOFF 

Plotter  Y-zero  position  offset 

0,100 

START 

Index  of  first  abscissa  value  to  plot 

1-10001 

UX1,UX2 

Absolute  Plotter  scaling  coordinates. 

0-16000 

UY1,UY2 

USER  size  plot  dimensions 

0-10100 

X1,X2 

Plotter  X  &  Y  coordinates  for  reading 

0-16000 

Y1,Y2 

soft  clip  windowing  points 

0-10100 

XP1,XP2 

Absolute  Plotter  scaling  coordinates. 

0-16000 

YP1,YP2 

for  WRITE  to  plotter 

0-10100 

NDATA(2) 

Number  of  data  points  in  spectrum 

1-10001 

PARAM(2,49) 

Instrument  operating  modes  table 

1-16 
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Table  III 


Glossary  of  REAL  Variables 


Name 


Description 


BAND 

BHEIGHT 

BWIDTH 

COFF 

ESCALE 

EMULT 

EOFF 

FXL.LXL 

FXT.LXT 

FYL.LYL 

FYT.LYT 

GAIN 

HEIGHT 

LHEIGHT 

LWIDTH 

LMIN 

LMAX 

NUMBER 

PMIN 

PMAX. 

RATIO* 

STEP 

T1 . T4 

WIDTH 

WN,WX 

WMIN,XMAX 

X,Y 

Xlabel 

XMIN.XMAX 

XMULT 

XOFF 

Xtick 

XN,XX 

YN.YX 

Ylabel 

YMULT 

YOFF 

YSCALE 

Ytick 


Spectral  Bandwidth  (nm)  -  AUTO  GAIN  mode 

Large  label  character  height 

Large  label  character  width 

Character  offset  variable,  labelling  mode 

Y-axis  scale  length.  Ext.  Coeff.  mode 

Power  of  10  for  scaling.  Ext.  Coeff.  mode 

Y-axis  offset.  Ext.  Coeff.  mode 

First  &  last  X-axis  labels 

First  &  last  X-axis  tick  marks 

First  &  last  Y-axis  labels 

First  &  last  Y-axis  tick  marks 

Instrument  gain  -  AUTO  SLIT  mode 

Normal  label  character  height 

Small  label  character  height 

Small  label  character  width 

Current  X-axis  lower  limit,  data  entry  check 

Current  X-axis  upper  limit,  data  entry  check 

General  purpose  data  entry  variable 

Pen  scale  minimum  limit 

Pen  scale  maximum  limit 

General  purpose  scaling  variable 

Step  size  (nm)  interval,  plotting  mode 

Powers  of  10  (10,...,1E4) 

Normal  label  character  width 

Current  plot  limits,  default  =  WMIN.WMAX 

Absolute  scan  limits  (nm)  or  (kK) 

General  purpose  plotting  variables 

X-axis  label  spacing 

Left  &  Right  edges  of  X-scale  on  plot 

Multiplier  to  scale  user  X-scale  to  10000  units 

Offset  to  user  X-scale  units  to  set  min'm  at  0 

X-axis  tick  spacing 

Data  entry  variables  for  X-scale,  ticks  &  labels 
Data  entry  variables  for  Y-scale,  ticks  &  labels 
Y-axis  label  spacing 

Scales  Y-axis  data  &  labels  for  Ext.  Coeff.  mode 
Offset  variable  on  current  Y-axis  scale 
Current  Y-axis  scale  length 
Y-axis  tick  spacing 


Table  IV 


Glossary  of  REAL  Array  Variables 


Name 

Description 

ORD(2) 

Final  ordinate  value  in  file 

ABSC(2) 

Final  abscissa  value  in  file 

CELL(2) 

Final  cell  #  value  in  file 

CYCLE(2) 

Final  cycle  #  value  in  file 

SAMPLE(2) 

Final  sample  #  value  in  file 

WAVE(2) 

Final  wavelength  value  in  file 

TIMER(2) 

Final  time  value  in  file 

DIST(2) 

Final  distance  value  in  file 

A(2. 10001) 

Absorbance  or  %T  array 

CONC(2) 

Concentration  of  sample  (M),  from  data  file 

FACTOR(2) 

Concentration  •  Pathlength  for  scaling  spectra 

FMAX(2) 

Starting  wavelength  of  scan  (nm) 

FMIN(2) 

Ending  wavelength  of  scan  (nm) 

FSTEP(2) 

Step  size  (nm),  from  data  file 

PATH(2) 

Pathlength  of  sample  cell  (cm),  from  data  file 

VARIABLE(2,14) 

Instrument  operating  conditions  table 

W(2,10001) 

Wavelength  array 
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Table  V 


Glossary  Of  CHARACTER  Variables 


Name 

1 

Description 

- 

1 

Screen  Control: 

1 

BELL 

CHAR(7)  bell  character 

CLR*2 

Clear  screen 

1  DOWN*2 

Move  cursor  down  1  line 

ESC 

CHAR(27)  escape  character 

HOME*2 

Move  cursor  to  upper  left  comer 

UP»2 

Move  cursor  up  1  line 

File  Status: 

ft 

DATE(2)*8 

Date  (mm/dd/yy) 

Fname(2)*20 

Filename,  Cl  convention 

LABEL(2)*72 

Descriptor  of  spectrum 

Pstr(23,16)*l4 

Table  of  parameter  setting  names 

1  Smin(2)*4 

Spectrum  ending  wavelength  (nm) 

1  Smax(2)*4 

Spectrum  starting  wavelength  (nm.) 

Sinc(2)*4 

Step  size  interval  (nm) 

Sstat(2)*10 

Spectrum  status  (VALID,  INVALID) 

1  Program  Control: 

C 

Literal  comma 

Code*2 

Main  menu  selection,  valid  until  reset 

Directory*  40 

User  directory  name.  Cl  convention 

Ext*14 

Extinction  Coefficient,  d  gitize  &  label  modes 

j  Filespec*63 

Full  pathname.  Cl  convention 

i 

Icode 

General  purpose  selection  key  entry 

Mult*5 

Y-axis  multiplier  label.  Ext.  Coeff.  mode 

Pcode*2 

Plotting  menu  option,  valid  until  reset 

Psize*8 

Screen  label  for  axes  type  selected 

Size 

Select  type  of  plotted  axes,  'N',  'F'  or  'U' 

^  String*!  4 

String  to  pass  data  to  or  from  subroutines 

Text*40 

General  purpose  string  for  labelling  plot 

TITLE*72 

String  to  be  printed  to  screen 

Xcode 

Selects  X-axis  scale,  (nm)  or  (kK) 

Xord*7 

X-axis  label  (nm)  or  (cm“*),  digitize  mode 

Ycode 

4 

Selects  Y-axis  scale,  'A',  ’E',  'R'  or  'T' 

» 

i 

4 
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...cont'd 


Plotter  Control: 

CR 

Etx 

LF 

Length*  6 

Lxnode*8 

Pattem*6 

Plotter*10 

Vaxes*2 

Vspec*2 

Xleft*6 

Xright*6 

Xstring*15 

Ylow*6 

Ystring*!  5 

Ytop*6 


CHAR(13),  carriage  return 

CHAR(3),  Terminator  for  label  mode  on  plotter 

CHAR(IO),  line  feed 

Broken  line  style  pattern  length  (%) 

Line  type  mode  (SOLID  or  BROKEN) 

Broken  line  type  number 

Uteral  'HP  9872A’  or  ’HP  7550A'  identifier 

Pen  velocity  (cm/ sec),  axes  vectors 

Pen  velocity  (cm/sec),  spectrum 

Literal  'O'  or  '-10000',  for  scale  command  SC 

Literal  '10000*  or  'O',  for  scale  command  SC 

X-axis  label,  plotter  output 

Literal  '0* 

Y-axis  label,  screen  and  plotter  output 
Uteral  '10000' 


1  FTN7X.L 

2  $FILES  0,1 

3  $EMA/DATA/ 

4  PROGRAM  PLOTS PEC 

5  c  - 

6  C 

7  C  This  Program  Is  Designed  To  Plot  Spectral  Data  Acquired  From 

8  C  The  CARY  2390  UV-VIS-NIR  Spectrophotometer  U/ith  CARYSPEC.RUN 

9  C 

10  C  The  Program  Supports  Output  To  HP  9872A  And  HP  7550A  Plotters 

11  C  The  HP  9872a  Is  Plotter  #1  At  Address  LU  33  (IEEE-488  Device  «2) 

12  C  The  HP  7550A  Is  Plotter  #2  At  Address  LU  17  (RS-232C  Device) 

13  C 

14  C  - 

15  C 


16 

C 

AUTHOR : 

17 

C 

18 

c 

19 

c 

20 

c 

21 

c 

WRITTEN: 

22 

c 

23 

c 

VERSION; 

24 

c 

25 

c 

REVISED: 

26 

c 

27 

c 

28 

c 

29 

c 

30 

c 

31 

c 

32 

c 

33 

c 

34 

c 

35 

c 

36 

c 

37 

c 

38 

c 

39 

c 

40 

c 

41 

c 

42 

c 

43 

c 

44 

c 

45 

c 

46 

c 

h7 

c 

48 

c 

49 

c 

50 

c 

C  1 

c 

52 

c 

53 

c 

Dr.  Robert  A.  Binstead, 
Chemistry  Division,  Code  6125, 
Naval  Research  Laboratory, 
Washington.  D.C.  20375-5000 


February,  1987 
1.72 

March,  1987: 


-  Debugged  plotting  EX)  loops 

-  Added  Ordinate/Abscissa  mode  checks 

-  Added  ending  wavelength  check 


Apr i 1 ,  1987 ; 

-  Altered  to  use  expc ' 1  X-values 

-  Added  screen  display  routine 

-  Removed  rounding  from  SUB  Xlimits 

-  Altered  rounding  to  use  NINT  function 


May,  1987; 


-  Altered  annotation  mode  to  support 
Title,  Parameters  &  Labels  on  plot 

-  Altered  AXES  selection  to  provide 
default  USER  values  &  easy  updates 

-  Changed  FFRCL  to  79  characters  to 
prevent  unwanted  line  wrapping 


June ,  1987 : 


Altered  Filename  convention  to  match 
the  directory  structure  of  the  new  Cl 
operating  system 

The  program  had  to  be  stripped  down 
to  run  under  the  Cl  operating  system 
by  eliminating  most  I/O  error  checks 


July,  1987: 


-  SEGMENTED  versions  of  the  full  size 


program  failed  with  Memory  Protect 
errors  at  the  INQUIRE  statement  if 
the  file  existed  already 


54  C 

55  C 

56  C 

57  C 

58  C 

59  C 

60  C 

61  C 

62  C 

63  C 

64  C 

65  C 

66  C 

67  C 

68  C 

69  C 

70  C 

71  C 

72  C 

73  C 

74  C 

75  C 

76  C 

77  C 

78  C 

79  C 

80  C 

81  C 

82  C 

83  C 

84  C 

85  C 

86  C 

87  C 

88  C 

89  C 

90  C 

91  C 

92  C 

93  C 

94  C 

95  C 

96  C 

97  C 

98  C 

99  C 

100  C 

101  C 

102  C 

103  C 

104  C 

105  C 

106  C 


August,  1987: 

-  Eliminated  INQUIRE  statement  for  file 
exists  or  open  check  to  avoid  system 
crash  in  the  segmented  version  of  the 
program.  The  error  checks  are  now  made 
using  error  numbers  returned  in  the 
OPEN  statement . 

Apr i 1 ,  1988 : 

-  Removed  needless  DATA  for  Pstr,  unused 
arrays  and  most  I/O  error  checks  on 
WRITE  to  allow  room  for  the  complete 
program  in  a  single  32K  word  segment. 

-  Added  Uppercase  string  conversion 
routine  for  all  string  entries. 

-  Subroutines  Xlimits  &  Ylimits  were 
replaced  by  Subroutine  Order (MIN , MAX) 

-  Subroutine  Xaxis  was  added  to  avoid 
duplication  in  setting  X-axis  Limits, 
Xtick  and  XIabel  spacing. 

-  Subroutine  Yaxis  was  added  to  avoid 
duplication  in  setting  Y-axis  Limits, 
Ytick  and  Ylabel  spacing. 

-  Altered  Annotation  Mode  For  Labels  Up 
To  40  Characters  In  Centered,  Left  & 
Right  Justified  Formats.  Positioning 
&  Character  Size  Commands  Were  Added. 

-  The  Pen  Position  Digitizing  Routine 
Was  Rewritten  In  Subroutine  Form  And 

A  Labelling  Command  Added  To  The  Menu. 

-  The  Annotation  Routine  Was  Augmented 
To  Support  Peak  Position  Labels  In 
Centered,  Left  &  Right  Justified  Modes 
With  Special  Offset  Features  To  Avoid 
Writing  Over  The  Spectrum. 

-  %  Transmission  &  %  Reflectance  Modes 
Were  Added  To  The  Y-scale  Selections. 

-  Re-scaling  Between  %T  &  Absorbance 
Was  Added  To  Support  Both  Plotting 
Modes  From  Either  Type  Of  Data  File. 

-  Rounding  Functions  Were  Altered  To 
Use  ANINT  (Whole  Number)  Command  To 
.‘.void  Integer  Overflow  With  ^'^T  &  WR 
Labe  1 1 i ng . 

-  Added  Automatic  Tick  &  Label  Routines 
To  Select  Preferred  Values  Based  On 

X  &  Y  Scale  Ranges.  Defaults  To  '.,10 
Of  Scale  For  Entry  Of  An  Anomalously 
Small  Or  Large  Range. 


-  Added  Code  For  Selection  Of  Plotter 
Type  (HP  9872A  /  HP  7550A)  And  Made 
Adjustments  To  The  P1,P2  Scaling 
Points  To  Fit  Output  Within  The  Hard 
Clip  Limits  Of  Both  Plotters. 

Y-AXIS;  Absorbance  or  Extinction  Coefficient 

X-AXIS;  Wavelength  or  Wavenumbers  x  lE-3 

MEMORY;  30,000  Words  (PROGRAM)  +  80,000  Words  LMA  (DATA) 


I NTEGER  D I FF , EX I T , MENU , PLOT , READ , S ETT I NGS , PARAM (2,49) 

INTEGER  FINISH, IX, lY, IZ,NCOL,ND,NDATA(2) ,NP , NV , PN , PU , PXOFF , PYOFF 
INTEGER  START,UX1 ,UX2,UY1 ,UY2,X1 ,X2 , Y1 , Y2 , XPl ,XP2,YP1 ,YP2 
REAL  BAND , COFF , ESCALE , EMULT , EOFF , CA I N ,  LM I N ,  LMAX , NUMBER 
REAL  PMIN,PMAX,RAT10,SC.'VLE,STEP,V,WN,WX,WMIN,WWAX 
REAL  XM I N , XMAX , XMULT , XOFF ,X,XN,XX,Y,YN,YX, YMULT , YOFF , YSCALE 
REAL  Xt ick,FXT,LXT,Yt ick,FYT,LYT,Xlabel , FXL , LXL , Y1 abe 1 ,FYL,LYL 
REAL  T1 , T2, T3, T4, LWIDTH, WIDTH, BW1DTH,LHEIGHT, HEIGHT, BHEICHT 
REAL  ABSC(2) ,CELL(2) ,CONC(2) ,CYCLE(2) ,D1ST(2) ,ORD(2) ,PATH(2) 

REAL  SAMPLE(2) ,TIMER(2) ,TEMP(2) ,WAVE(2) ,FMIN(2) ,FMAX(2) ,FSTEP(2) 
REAL  FACTOR(2)  ,.VARI ABLE(2 , 14) ,  A(2 , 10001 ),  W(2 , 10001 ) 

LOGICAL  MATCH 

Dimension  Screen  Control  String  Variables 
CHARACTER  BELL , CLR*2 , DOWN*2 , ERASE*2 , ESC , HOME*2 , UP* 2 

Dimension  Program  Parameter  Variables 

CHARACTER  DATE(2) *8 , Fname (2)*20 , LABEL( 2)*72 , PI ot t er ( 2) *1 0 
CHARACTER  Smi n(2)*4 , Smax(2)*4 , S i nc (2)*  ,Sstat(2)*10 
CHARACTER  C, CR, Code*2 . Di rectory*40 ,Etx,Ext*I4,Fi I espec*63 
CHARACTER  I  code , Lengt  h*6 , LF, Lmode*8 , Mul t*5 , Pcode*2 , Pat  t  ern*6 
CHARACTER  Ps i ze*8 , S i ze , St  ring* 14 , Text *40 , TITLE*? 2 , Vaxes*2 
CHARACTER  Vspec*2 , Xcode ,Xle ft*6 .Xord*7 ,Xright*6,Xstring*15 
CHARACTER  Ycode , Yl ow*6 , Yst r i ng*l 5 , Yt op*6 , Ps t r ( 23 , 16)*14 


COMMON  /DATA/A, W 

COMMON  /XPARAM/LMIN.LMAX.WN.WX.Xt ick.FXT.LXT.Xlabel  , FXL, LXL 
COMMON  /YPARAM/Yt ick,FYT,LYT,Ylabel .FYL.LYL 

COMMON  /DIG IT/FACTOR, XMIN, XMAX, YOFF, YSCALZ, IX, lY, J ,PU,X1 ,X2,Y1  ,Y 


Define  Rounding  Functions  Using  Nearest  Integer  Commands 


R0(X)-N1NT(X) 


!  Tvpe  =-  INTEGER 


160  R1(X)-ANINT(X*10.0)/10.0  !  Type  -  REAL 

161  R2(X)-ANINT(X*100.0)/100.0  !  Type  -  REAL 

162  R3(X)-ANINT(X*1.0E3)/1.0E3  !  Type  -  REAL 

163  C 

164  C  - 

165  C 

166  C  Initialize  String  Variables 

167  C 

168  DATA  (Plotter(I) , I-l .2)/'HP  9872a  ','HP7550a  '/ 

169  DATA  (Pstr(l , I) , I-l ,6)/'ABSORBANCE’ TRANSMISSION' , 

170  & 'TEMPERATURE' , REFLECTANCE' , 'CONCENTRATION' ,' EMISSION' / 

171  DATA  (Pstr(2. I) . I-l ,4)/'WAVELENGTH' , 'TIME' , 'TEMPERATURE' , 

172  &' DISTANCE'/ 

173  DATA  (Pst r (3, I). 1-1, 11)/' OFF’ , ’ 0 . 01 ' , ' 0 . 02 '  , ' 0 . 05 '  , ' 0 . 1'  ,  ' 0 . 2 '  , 

174  &'0.5' , ' 1 .0' , '2.0' , '5.0' , ' 10.0'/ 

175  DATA  Pstr(4, l)/'OFF'/ 

176  DATA  (Pstr(4. I) , I-6,15)/’0.2’ , '0.5' . ' 1 .0' , '2.0' , '5.0' , 

177  &’ 10’ . '20' , '50' . ’ 100’ , '200'/ 

178  DATA  (Pstr(5, I) , I-l, 4)/'AUT0  SELECT' , 'AUTO  CAIN' , 'AUTO  SLIT' , 

179  &' SINGLE  BEAM'/ 

180  DATA  (Pstr(6. I) , 1-1 ,5)/'OFF’ . 'NORMAL' , ’ 1ST  DER1V','2ND  DERIV  , 

181  &'LOG'/ 

182  DATA  (Pstr(7,I),l-l,9)/'0.01' .'0.02' , '0.05'  ,  'O.r  ,  '0.2'  ,  '0.5'  , 

183  &' 1 .0' , '2.0' , '4.0'/ 

184  DATA  (Pstr(8, I) , I-IO, 16)/’2’ . '5' , ' 10' , '20' , '50' , ' 100' , ' 200'/ 

185  DATA  (Pstr(9,I),I-12,15)/’10','20','50' ,'100'/ 

186  DATA  (Pstr(10,I),l-l,5)/’-1.9  TO  0.6’, '-2.0  TO  0.5’, 

187  &'-2.1  TO  0.4' ,’ -2.2  TO  0.3' ,’ -2.3  TO  0.2'/ 

188  DATA  (Pstrdl  ,  I)  ,  I-l  ,  16)/' +/-0 . 01  ’  ,  ’+/-0.02’  ,  ’+/-0.05’-,  ’+/-0. 1  '  , 

189  &'+/-0.2' , '+/-0.5' , '+/-1-0' , '+/-2.0' , '+/-5.0' , '+/-10' . ’+/-20’ , 

190  &'+/-50' , '+/-100' , '+/-200' . ’+/-500’ , ’+/-1000'/ 

191  DATA  (Pstr(12, I) , I-l , I6)/'+/-0 . 01 ' , '+/-0.02' , '+/-0.05' , '+/-0. i ' , 

192  &'+/-0.2' , '+/-0.5' , '+/-^ -O' , '+/-2.0’ , ’+/-5.0' , '+/-10' , ’+/-20’ , 

193  &'+/-50' , '+/-100' , ’+/-200’ , '+/-500' , '+/-1000’/ 

194  DATA  (Pstr(14,  I)  ,  I-l ,  10)/'0' ,  ’  10’  '  20 '  ,  '  30 '  ,  '  40 '  ,  '  50 '  ,  '  ..0  ’  ,  '  70 '  , 

195  &'80','90'/ 

196  DATA  (Pstr(15, 1) , I-l ,4)/'0.5' , ' 1 .0’ , ’ 3.0' , ' 10’ / 

197  DATA  (Pst r (16, I) , I-l ,2)/' NORMAL’ , 'REVERSE'/ 

198  DATA  (Pstr(17, I) , I-l ,2)/'OFF’ , 'ON'/ 

199  DATA  (Pstr(18, I) , I-l .2)/'REPEAT  SCAN' , 'SCL/MULTI '/ 

200  DATA  (Pstr(19, I) , 1-1 ,2)/’SERIAL' , 'OVERLAY'/ 

201  DATA  (Pstr(20, I) , I-l ,4)/'BOTH  ON' , 'UV  ONLY' . 'VIS/NIR  ONLY' , 

202  &’BOTH  OFF’/ 

203  DATA  (Pstr(21 , I) , I-l ,3)/'AUTO' , 'UV' . 'VIS/NIR' / 

204  DATA  (Pstr(22, 1) , 1-1 ,3)/' AUTO' , ' UV/VIS’ , 'NIR' / 

205  DATA  (Pstr(23, I) , I-l ,2)/'FULL' , ' 1/3'/ 

206  C 

207  BELL-CHAR(7) 

208  ESC-CHAR(27) 

209  CLR-ESC//'J' 

210  HOME-ESC//'h' 

211  UP-ESC//'A' 

212  DOWN-ESC//' B' 

213  ERASE-ESC//’K' 


!  Select  Plotter  =»1  (HP  9872a) 


214  C 

215 

216 

217 

218 

219 

220  C 

221  C 

222  C 

223  C 

224  C 

225  C 

226  C 

227 

228 

229 

230 

231 

232 

233  C 

234  C 

235  C 

236  C 

237  C 

238  C 

239  C 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 


PN-1 

Tl-10.0 

T2-100.0 

T3-1000.0 

T4-10000.0 


Assign  Statement  Labels 


ASSIGN  10  TO  MENU 
ASSIGN  100  TO  READ 
ASSIGN  300  TO  SETTINGS 
ASSIGN  370  TO  PLOT 
ASSIGN  390  TO  DIFF 
ASSIGN  9000  TO  EXIT 


Data  Transfer  and  Plot  Control  Menu 


CALL  FFRCL(79)  !  Set  Fields  to  79  Characters 

10  WRITE  (1,*)  HOME,CLR,'_' 

NCOL-70 

TITLE-'Cary  2390  Spectral  Data  Plotting' 

CALL  Center(TITLE) 

WRITE  (1 . ' (T61 ,A2.A9) • )  UP. 'Rev:  1.72’ 

CALL  Line(NCOL) 

WRITE  (1,20)  'CODE' , 'FUNCTION' . 'STATUS’ , 'MIN’ , 'MAX' INC 
"’3  FORMAT  (T4,A4,T18,A8,T38.A7,T50.A3,T58.A3.T66,A3) 

CALL  Line(NCOL) 

WRITE  (1,30)  'Rl',' . Read  Spect  rum  . ',Sstat(l), 

&Smi n( 1 ) , Smax( 1 ) , S i nc( 1 ) 

WRITE  (1,30)  'R2’  ,' . Read  Spectrum  *2 . '  ,Sstat(2)  , 

&Smi n(2) , Smax(2) , S i nc(2 ) 

WRITE  (1,50)  'IS',' . Instrument  Settings . ' 

WRITE  (1,30)  ’PS’.' . Plotter  Selection . '  ,  PI  ot  t  er  ( PN)) 

WRITE  (1,40)  'PI',' . Plot  Spectrum  =*1 . Fname  ( 1 ) 

WRITE  (1,40)  'P2'.’ . Plot  Spectrum  =2 . ',Fname(2) 

WRITE  (1,50)  '  DS  '  ,  ' . Difference  Spectrum . ' 

WRITE  (1,50)  'EX',' . EXIT  Menu . ' 

30  FORMAT  (/ ,  T5  ,  A2  ,  T7  .  A30  .  T39  ,  Al  0  ,  T50  ,  A4  ,  T58  ,  A4  ,  T56  ,  A4;i 
40  FORMAT  (/ ,  T5  ,  A2  ,T7  ,  A30  ,  T39  .  Al6'i 
50  FORMAT  (/ , T5 , A2 , T7 , A30 ) 

WRITE  (1,*) 

CALL  Line(NCOL) 

WRITE  (1,*) 

70  WRITE  (1,*)  UP, ERASE, 
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267 

268 

269 

270 

271  C 

272 

273 

274 

275 

276 
111 

278 

279 

280 
281 
282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297  C 

298  C 

299  C 

300  C 

301  C 

302  C 

303  C 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 


WRITE  (1 . ' (T3.A15,A.A2) • )  'Enter  the  CODE;' 
READ  (1,80)  Code 
80  FORMAT  (A2) 

CALL  Upper (Code) 


,BELL, 


IF  ((Code. EQ. 'R1 ') .OR. (Code. EQ. ■R2' ) )  THEN 
J-ICHAR(Code(2:2))-48 
GO  TO  READ 
END  IF 

IF  (Code.EQ. ' IS' )  GO  TO  SETTINGS 
IF  (Code.EQ. ' PS ' )  THEN 

WRITE  (1,*)  UP, ERASE,'  Select:  1 . . . HP  9872A,' 
it  '  2.  .  .HP  7550A  ?  '  ,BELL,  '_' 

READ  (1,99)  I  code 

IF  ((Icode.NE. ' 1 ' ) .AND. (Icode.NE. '2' ) )  THEN 
CO  TO  90 
END  IF 

PN- 1 CHAR ( I  code) -48 
GO  TO  MENU 
END  IF 

IF  ((Code.EQ. 'PI' ) .OR. (Code.EQ. 'P2’ ))  THEN 
J-ICHAR(Code(2;2))-48 
GO  TO  PLOT 
END  IF 

IF  (Code.EQ. 'DS' )  THEN 
J-1 

CO  TO  DIFF 
END  IF 

IF  (Code.EQ. 'EX' )  GO  TO  EXIT 
CO  TO  70 


Read  Spectrum'  (J  -  1  or  2) 


99  FORMAT  (Al) 

100  IF  (Sstat(J) .EQ. 'VALID' )  THEN 

WRITE  (1,*)  UP, ERASE,'  Spectrum  is  PRESENT; 
&  'Proceed  (Y  or  N)  ?  ',BELL,'_' 

READ  (1,99)  I  code 
CALL  Uppe  r ( I  code ) 

IF  (Icode.EQ. 'N' )  GO  TO  70 
IF  (Icode.NE. 'Y' )  CO  TO  100 
END  IF 

WRITE  (1,*)  HOME,CLR,'_' 

TITLE-' Read  Spectrum  //Code (2 ; 2) 

CALL  Center(TITL£) 

CALL  Line(NCOL) 

WRITE  (1,*)  DOWN,'  Enter  Filename;  ',BELL.'_' 
READ  (1 , ' (A20) ' )  Fname(J) 

CALL  Upper(Fname(J) ) 
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320  WRITE  (1.*)  DOWN,’  Directory,  (Return  -  /DEFAULT/)  ’.BELL, 

321  READ  (1,'(A40)')  Directory 

322  IF  (Directory. EQ. '  ')  THEN 

323  Fi lespec-Fname( J) 

324  CO  TO  no 

325  END  IF 

326  K-40 

327  DO  WHILE  (Di rectory(K:K) .EQ. '  ') 

328  K-K-l 

329  END  DO 

330  IF  (Directory(K:K) .EQ. '/' )  R-K-l 

331  Fi lespec-Di rectory( 1 ; K) //' /' //Fname ( J ) 

332  no  K-63 

333  DO  WHILE  (Fi lespec (K: K) . EQ. ’  ') 

334  K-K-1 

335  END  DO 

336  WRITE  (1,*)  UP, ERASE,'  Validating;  ' , F i 1 espec ( 1 : K) , ' 

337  CALL  Wait (1.0) 

338  OPEN  (UNIT-66, FILE-Fi lespec(l ;K) , IOSTAT-N,STATUS-'OLD' ) 

339  WRITE  (1,*) 

340  IF  (N.NE.O)  THEN 

341  WRITE  (1,*)  UP, ERASE, 

342  N-N-500 

343  IF  (N.EQ.6)  WRITE  (1,*)  '  File  does  NOT  EXIST:  ’,BELL,'_' 

344  IF  (N.EQ.8)  WRITE  (1,*)  ’  File  is  already  OPEN:  ',BELL,'_ 

345  IF  ((N.NE.6)  .AND.  (N.NE.8))  WRITE  (1,*)  ’  Disk  Error  :»  ',N 

346  &  BELL,'  _' 

347  CLOSE  ( UN IT-66,  STATUS-’ DELETE ' ) 

348  CALL  Wait (2.0) 

349  WRITE  (1 ,*) 

350  Sstat (J)-' INVALID’ 

351  CO  TO  MENU 

352  END  IF 

353  WRITE  (1,*)  UP, ERASE,’  Reading  File;  F i 1 espec ( 1 : K) 

"34  READ  (66,FMT-120,  lOSTA"''  N,ERR-210)  LABEL(J) 

355  120  FORMAT  (A72) 

356  WRITE  (1,*)  DOWN,'  Title:’ 

357  WRITE  (1,*)  ’  ’ ,LABEL(J) 

358  READ  (66,FMT-130)  DATE(J) 

359  130  FORMAT  (A8) 

360  WRITE  (1,*)  DOWN,'  Date:  ’,DATE(J) 

361  READ  (66,*)  FMIN( J) , FMAX( J) , FSTEP( J ) . 

362  6cCONC(J)  ,PATH(J) 

363  READ  (66,*)  ORD(J) ,ABSC(J) ,CELL(J) ,CYCLE( J) , 

364  &SAMPLE(J) ,WAVE(J) ,TIMER(J) ,TEMP(J) ,D1ST(J) 

365  WRITE  (1,140)  DOWN ,’ Scan  Range :  ’,F\iAX(J),’  to  ’,FMIN(Jn 

36^  &’  nm  at  ’,FSTEP(J),’  nm  steps' 

367  WRITE  (1 , ’ (T4,A2,A12,2X,G9.4) ’ )  DOWN, 'Cone.  (M)  ;  ',CONC(J' 

368  WRITE  (1,150)  DOWN, 'Path  (cm)  ;  ’,PATH(J) 

369  WRITE  (1,150)  DOWN, 'Temp.  (C)  :  ',TEMP(J) 

370  140  FORMAT  (T4 , A2 , A12 , F6 . 2 , A4 , F6 . 2 , A7 , F4 . 2 , A9) 

371  150  FORMAT  (T4 , A2 , Al 2 , F6 . 2) 
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372 

373 

374 

375 

376 

377 

378 

379 

160 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

170 

391 

392 

393 

394 

395 

180 

396 

397 

190 

398 

399 

400 

401 

402 

403  C 

404 

200 

405 

406 

407 

210 

408 

409 

220 

410 

411 

412  C 

413 

230 

414 

415 

416 

417 

418 

240 

419 

420 

421 

250 

422 

423 

CALL  Str(FMIN(J) .String, 4) 

Smin(J)-String(2;5) 

CALL  S t r ( FMAX( J ) , S t r i ng , 4) 

Smax(J)-St  r ing(2 ; 5) 

CALL  Str(FSTEP(J) .String, 4) 

Si nc(J) -String (2:5) 

READ  (66,FMT-160. IOSTAT-N,ERR-210)  NP , NV , NDATA( J ) 

FORMAT  (13,13, 16) 

IF  ((NP.NE.49) .OR. (NV.NE. 14))  CO  TO  200 
IF  (ABS(FMIN(J)-WAVE(J)) .CT.O. 1)  THEN 

WRITE  (1,*)  DOWN,'  Warning!  Scan  ended  at ' , WAVE( J ) , '  nm' 
CALL  Wait (2.0) 

END  IF 
WRITE  (1,*) 

TITLE-'#^  READING  SPECTRUM  ###' 

WRITE  (1 ,*)  HOME, 

CALL  Center(TITLE) 

READ  (66,FMT-170)  (PARAM( J , K) , R-1 , NP) 

FORMAT  (12) 

READ  (66,*)  (VARIABLE(J .K) ,K-1 ,NV) 

READ  (66,*)  (A(J,k:),R-l,NDATA(J)) 

READ  (66,*)  (W(J ,K) ,K-1 ,NDATA(J) ) 

Sstat (J)-' VALID' 

CLOSE  (UNIT-66, IOSTAT-N,ERR-210) 

IF  (Sstat (J) .NE. 'VALID* )  GO  TO  MENU 

WRITE  (1,*)  UP, ERASE,'  EDIT  Cone. /Path:  (Y  or  N)  ?  '.BELL, 
READ  (1,99)  I  code 
CALL  Uppe  r ( I  code ) 

IF  (Icode.EQ.'Y')  GO  TO  230 
IF  (Icode.NE. 'N')  GO  TO  190 
GO  TO  MENU 


WRITE  (1,*)  UP, ERASE,'  Data  Format  Error:  ',BELL,'_' 
Sstat (J)-' ABORTED' 

CO  TO  220 

WRITE  (1,*)  UP, ERASE,'  Disk  Error  =*',N,BELL,' 

Sstat (J)-' INVALID' 

CALL  Wait (2.0) 

WRITE  (1,*) 

GO  TO  180 


WRITE  (1,*)  UP, ERASE,'  (C...Conc.  ,  P...Path)  ?  ',BELL,'_' 
READ  (1,99)  I  code 
CALL  Upper (1  code) 

IF  (Icode.EQ. 'P' )  CO  TO  250 
IF  (Icode.NE. 'C )  CO  TO  230 

WRITE  (\*)  UP, ERASE,'  Concentration,  (.M)  :  ',BELL,'_' 

READ  (1 ,*,ERR-240)  CONC(J^ 

CO  TO  190 

WRITE  (1,*)  UP, ERASE,'  Pathlength,  ic.m);  ',BELL,'_' 

READ  (1 ,*,ERR-250)  PATH(J) 

CO  TO  190 


.BELL 


•J 


424  C 

425  C  - 

426  C 

427  C  Error  Messages  &  Returns  To  Menu 

428  C 

429  C  - 

430  C 

431  260  TITLE-'  Cone,  or  Path  -  0.0  ' 

432  CO  TO  290 

433  270  TITLE-'  Spectra  are  MISMATCHED:  ' 

434  GO  TO  290 

435  280  TITLE-'  Spectrum  is  ABSENT:  ' 

436  290  WRITE  (1,*)  UP . ERASE ,TITLE( 1 : 26) , BELL, ' _' 

437  CALL  Wait (2.0) 

438  WRITE  (1 ,*) 

439  CO  TO  70 

440  C 

441  C  - 

442  C 

443  C  Display  Instrument  Settings 

444  C 

445  C  - 

446  C 

447  300  IF  ((Sstat(l).NE. 'VALID'). AND. (Sstat(2).NE. 'VALID'))  CO  TO  280 

448  310  WRITE  (1,*)  UP, ERASE,'  Spectrum  #,  (1  or  2):  ',BELL,'_' 

449  READ  (1 , ' ( 12) ' , ERR-310)  J 

450  IF  ((J.NE. 1) .AND. (J.NE.2))  CO  TO  310 

451  IF  (Sstat(J) .NE. 'VALID' )  GO  TO  280 

452  BAND-VARIABLE(J,10) 

453  GAIN-VARIABLE(J ,6) 

454  St r i ng-Pst r (7 , PARAM( J , 7)+l ) 

455  K-1 

456  IF  (PARAM(J, 1) .NE.O)  THEN 

457  String-Pstr(8,PARAM(J,8)+l) 

458  K-11 

459  END  IF 

460  CALL  Val (String, PMAX) 

461  PMIN-VAR1ABLE(J ,K) 

462  PMAX-PMIN+PMAX 

463  K-11 

464  IF  (PARAM(J ,6) .EQ.4)  K-10 

465  String-Pstr(K,PARAM(J,K)+l) 

466  WRITE  (1,*)  HOME,CLR,'_' 

467  TITLE— ' Instrument  Settings' 

468  CALL  Center(TITLE) 

469  CALL  Line (50) 

470  WRITE  (1  ,  '  (T20,A8,T40,A'')  '  )  '  FUNCTION'  ,' SETTING' 


471  CALL  LINE(50) 

472  WRITE  (1 ,*) 

473  WRITE  (1,320)  'WAVELENGTH  LIMITS . FMIN(J FMAX(J ) 

474  WRITE  (1,330)  'ORDINATE . '  ,  Pst  r  ( 1  ,  PARA'^(  J  ,  1  ) +1 ) 

475  WRITE  (1,330)  'ABSCISSA . '  ,  Ps  t  r  ( 2  ,  PARA^C  J  ,  2 ) -rl ) 

476  WRITE  (1,330)  ' SCAN  RATE  (nm/sec) . ' , Ps t r ( 3 , PARAM( J , 3 ) +1) 
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477 

478 

479 

480 

481 

482 

483 

484 

485 

486 

487 

488 

489 

490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500  C 

501  C 

502  C 

503  C 

504  C 

505  C 

506  C 

507 

508 

509 

510 

511 

512 

513 

514 

515 

516  C 

517 

518 

519 

520 

521 

522 

523 

524 

525 

526 

527 

528 

529 


WRITE  (1,330)  'CHART  DISPLAY  (nm/cm) . . ' , Pst r (4 , PARAM( J , 4)+l) 

WRITE  (1,330)  'REFERENCE  MODE . '  ,  Pst  r  (5  ,  PARAM(  J  ,  5 ) +  1 ) 

WRITE  (1.340)  'SBW  (nm)  .  CAIN . '  ,  BAND,  '  ,  '  ,  C.  .IN 

WRITE  (1,330)  'PEN  FUNCTION . '  ,  Pst  r  (6  ,  PARAM(  J  .  6)  +  l ) 

WRITE  (1,340)  'PEN  LIMITS  (Min, Max) . . . ' , PMIN , ' , ' , PMAX 
IF  (PARAM(J,6) .CT. 1)  WRITE  (1,350)  UP, ERASE, S t r i ng 

WRITE  (1,330)  'RESPONSE  TIME  (sec) _ ' , Pst r ( 15 , PARAM( J , 13)+1 ) 

WRITE  (1,330)  'BEAM  INTERCHANGE . ' , Ps t r ( 16 , PARAM( J , 16) +1 ) 

WRITE  (1,330)  'SLIT  HEIGHT . '  ,  Pst  r  (23  ,  PARAM(  J  ,  23)  +  l ) 

WRITE  (1,330)  'LAMP  SELECT . ’  ,  Pst  r  (2 1 .  PARAM(  J  .  2 1 )  + 1 ) 

WRITE  (1,330)  'DETECTOR  SELECT . '  , Pst r (22 . PARAM( J , 22) +1 ) 

320  FORMAT  (T15 , A23 ,T40 , F5 . 2 , A. F5 . 2) 

330  FORMAT  (T15 , A23 ,T40 ,A14) 

340  FORMAT  (T15 , A23 .T40 . F4 . 2 , A, F5 . 2) 

350  FORMAT  (T40 , A2 ,A2 ,A14) 

WRITE  (1,*) 

CALL  Line(50) 

WRITE  (1,*) 

360  WRITE  (1, ' (T14,A2,A2,A23,A,A)')  UP, ERASE, 

&'Press  RETURN  for  MENU;  ',BELL,'_' 

READ  (1,99)  I code 
IF  (Icode.NE. '  ')  CO  TO  360 
GO  TO  MENU 


WRITE  (1,330) 
WRITE  (1,330) 
WRITE  (1,330) 
WRITE  (1,330) 


Plot  Mode  Setup 


370  IF  (Sstat(J) .NE. 'VALID' )  GO  TO  280 
GO  TO  400 

390  IF  ((Sstat(l) .NE. 'VALID' ) .OR. (Sstat(2) .NE. 'VALID' ))  CO  TO  280 
MATCH- . TRUE . 

IF  (FSTEP(l) .NE.FSTEP(2))  MATCH-. FAL'o. 

IF  (FMIN(l) .NE.FMIN(2))  MATCH-. FALSE. 

IF  (FMAX(l) .NE.FMAX(2))  MATCH- . FALSE . 

IF  (MATCH)  CO  TO  400 
GO  TO  270 

400  TITLE-'Plot  Spectrum  =*' //CHAR(  J+48) 

FACTOR ( 1 ) -CONC ( 1 ) *P ATH ( 1 ) 

FACTOR ( 2 ) -CONC ( 2 ) *PATH ( 2 ) 

IF  (Code.EQ. 'DS' )  THEN 

IF  ((FACTOR(l) .EQ.0.0) .OR. (FACTOR(2) .EQ.0.0) )  CO  TO  260 
RATIO-FACTOR(l)/FACTOR(2)  !  Temporary  Scaling  Factor 

TITLE-' Di fference  Spectrum  Plot' 

END  IF 

IF  (FACTOR(J) .EQ.0.0)  FACTOR( J)-l . 0  !  Avoid  Division  By  Zero 
WRITE  (1 ,*)  HOME.CLR, 

CALL  Center(TITLE) 

CALL  Line(NCOL) 

WRITE  (1 ,*)  DOWN, DOWN 
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I 


530 

531 

532 

533 

534  410 

535  420 

536 

537 

538 

539 

540 

541 

542 

543 

544 

545 

546 

547 

548 

549 

550 

551 

552  430 

553 

554 

555 

556 

557  C 

558 

559 

560 

561 

562 

563 

564 

565 

566 

567  440 

568 

569 

570 

571  C 

572 

573 

574 

575  C 

576 

577 

578 

579 

580 

581 

582 

583 


WRITE  ( 1 , ' (TIO , A32) ' )  'X-scale  units;  N . Nanometres' 

WRITE  (1,*) 

WRITE  ( 1 , ' (T26 , A17) ' )  'W . Wavenumbers' 

WRITE  (1,*)  DOWN 

WRITE  (1,420)  UP. ERASE, 'Code;  '.BELL,'_' 

FORiMAT  (T19,A2,A2,A7,A,A) 

READ  (1,99)  Xcode 
CALL  Upper(Xcode) 

IF  (Xcode. EQ. 'W' )  THEN 
WRITE  (1,*)  UP. ERASE, UP, UP, ERASE, UP, UP 
WRITE  (1 . ' (T26.A2,A2,A11) • )  UP, ERASE, 'Wavenumbers ’ 

Xord-'  (cm-1)' 

WMAX-R2(T4/FMIN(J)) 

WMIN-R2(T4/FMAX(J) ) 

GO  TO  430 
END  IF 

IF  (Xcode. NE. 'N' )  CO  TO  410 
WRITE  (1,*)  UP, ERASE, UP, UP, ERASE, UP, UP 
WRITE  (1 , ' (T26.A2,A2.A10) • )  UP, ERASE, ’ Nanomet res ' 

Xord-'  (nm)' 

WMAX-R1(FMAX(J)) 

WMIN-R1(FMIN(J)) 

NI^NDATA(J)  !  Number  Of  Data  Points 

STEP— FSTEP( J)  !  Step  Size  (nm) 

XMIN-WMIN  !  X-scale  MIN/MAX  -  Scan  Limits 

XMAX-WMAX 

CALL  Xaxis(XMIN,XMAX)  !  Set  Limits,  Tick  &  Label  Spacing 


WRITE  (1,*)  DOWN, DOWN 

WRITE  ( 1 , ' (TIO , A32) ' )  'Y-scale  units;  A . Absorbance' 

WRITE  (1,*) 

WRITE  (1 , ' (T26,A28) ' )  'E . Ext i net i on  Coef f i c i ent ' 

WRITE  (1,*) 

WRITE  (1 , ' (T26.A19) ' )  'R . %  Reflectance' 

WRITE  (1,*) 

WRITE  (1 , ' (T26,A20) ’ )  'T . X  Transmission' 

WRITE  (1 ,*)  DOWN 

WRITE  (1,420)  UP, ERASE, 'Code;  ',BELL,'_' 

READ  (1,99)  Ycode 
CALL  Upper (Ycode) 

Text-UP//ERASE//UP//UP//ERASE//UP//UP//ERASE//UP//UP//ERASE//'UP 
IF  (Code.EQ. 'DS' )  THEN 

IF  ((Ycode. EQ. 'R' ) .OR. (Ycode .EQ. 'T' ) )  CO  TO  440 
END  IF 

IF  (Ycode. EQ. 'A' )  THEN 

Ystring-' ABSORBANCE'  !  Select  Y-axis  Label  For  Plot 

K-0  !  Ordinate  Mode  =  Absorbance 

ELSE  IF  (Ycode.EQ. 'E' )  THEN 

K-0  !  Ordinate  Mode  -  Absorbance 

ELSE  IF  (Ycode.EQ. 'R' )  THEN 
Ystring-'%  Reflectance' 

K-1  !  Qrdinate  Mode  -  ^  Transmission 
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584 

585 

586 

587 

588 

589 

590  C 

591 

592 

593 

594 

595 

596 

597 

598 

599 

600 
601 

602  C 

603 

604 

605 

606 

607 

608 

609 

610 
611 
612 

613  C 

614 

615 

616 

617 

618 

619 

620 
621 
622 

623  C 

624 

625 

626 

627 

628 

629 

630 

631 

632 

633 

634 

635 

636 


ELSE  IF  (Ycode.EQ. 'T' )  THEN 
Ystring-'%  Transmission' 

K"1  !  Ordinate  Mode  -  %  Transmission 

ELSE 

CO  TO  440 
END  IF 


IF  (PARAM(J , 1) .NE.K)  THEN  ! 

CALL  Rescale(J ,K,ND)  ! 

PARAM(J,1)-K  ! 

END  IF 

IF  (Code.EQ. 'DS' )  THEN  ! 

K-0  ! 

IF  (PARAM(2, 1) .NE.K)  THEN  ! 

CALL  Rescale(2.K.ND)  ! 

PARAM(2,1)-K  ! 

END  IF 
END  IF 


Test  Ordinate  Mode  Of  Spectrum  =J 
Rescale  Between  Absorbance  &  ‘M 
Reset  Ordinate  Variable  To  Match 

Difference  Spectrum  Mode 
Only  Absorbance  Mode  Allowed 
Test  Ordinate  .Mode  Of  Spectrum  =2 
Rescale  To  Absorbance  If  Required 
Reset  Ordinate  Variable  To  Match 


IF  (Ycode.NE. 'E' )  THEN 

WRITE  (1,*)  Text ( 1 : 24) , UP  !  Erase  Unwanted  Labels 

WRITE  (1 , ' (T26.A2.A2,A15) ' )  UP, ERASE ,Yst r i ng 

YMULT-l.O  !  YMULT  Scales  Ext.Coeff.  Mode 

IF  (Ycode.EQ. 'A' )  CO  TO  450 

YOFF-0 . 0 

YSCALE- 100.0 

CALL  Yax is (YSCALE, Ycode) 

CO  TO  460 
END  IF 


IF  (FACTOR(J) .EQ.0.0)  THEN 
WRITE  (1 , ' (T19,A2,A2,T30,A19,A,A) ' )  UP, ERASE, 

&  'Cone,  or  Path  —  0 . 0 ' , BELL, ' _' 

CALL  Wait (2.0) 

WRITE  (1,*) 

CO  TO  440 
END  IF 

WRITE  (1,*)  Text ( 1 ; 24) , UP  !  Erase  Unwanted  Labels 

WRITE  (1 , ' (T26,A2,A2,A22) ' )  UP, ERASE, ' Ext i net i on  Coe f f i c i ent 


450  K-1 

PMAX-0 . 0 

DO  WHILE  (K.LE.ND)  !  Find  Maximum  Absorbance 

NUMBER-A(J ,K) 

IF  (Code.EQ. 'DS')  THEN 

NUMBER-ABS(A(1 ,K) -A(2 , K)*RATIO) 

END  IF 

IF  (PMAX.LT. NUMBER)  PMAX-NUMBER 
K-K+1 
END  DO 

SCALE— R3  (P,MAX+0 . 001 )  !  Round  Up  Absorbance  Max'm 

IF  (PMAX.CT.0.01)  SCALE-R2(PMAX+0.01) 

IF  (PMAX.CT.O. 10)  SCALE-Rl (PMAX+0. 10) 
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Initial  Y-scale  Min'm 
Initial  Y-scale  Max'm 
Initial  E-scale  Max'm 


637  YOFF-0 . 0 

638  YSCALE-SCALE 

639  ESCALE-SCALE 

640  IF  (Code.EQ. 'DS' )  THEN 

641  YSCALE-2 . 0*SCALE  !  2  x  Y-scale  for  +/-  Max'm  Range 

642  YOFF-SCALE  !  Set  Zero  Position  At  Half  Y-scale 

643  END  IF 

644  CALL  Yaxis(YSCALE,Ycode)  !  Set  Tick  &  Label  Spacing 

645  IF  (Ycode.EQ. 'A' )  CO  TO  460 

646  C 

647  C  - - - 

648  C 


649  C 

Scale  YMULT  for  Extinction  Coefficient  Mode 

650  C 

651  C 

DO  WHILE  (ESCALE/FACTOR(J) 

.LT. 

1.0) 

652  C 

653 

654 

ESCALE-ESCALE*10 . 0 

655 

END  DO 

656 

DO  WHILE  (ESCALE/FACTOR(J) 

.CT. 

10.0) 

657 

ESCALE-ESCALE/ 10.0 

658 

END  DO 

659 

RAT 1 0-ESCALE/SCALE 

1 

Power  Of  10  For  Exponent 

660 

YMULT-RAT I O/FACTOR ( J ) 

1 

Scaling  Factor  For  Ext.  Coeff. 

661 

662  C 

663  C 

CALL  Exponent (RATIO, Mul t ) 

Convert  Exponent  To  String 

664  C 

665  C 

Plotting  Routine:  HP  9872A  &  HP  7550A 

666  C 

667  C 

668  C 

669  460  PU-33 

t 

Plotter  Logical  Unit  (HP 

9872A) 

670 

PXOFF-0 

1 

X-zero  Position  Offset  (HP 

9872A) 

671 

PYOFF-0 

r 

Y-zero  Position  Offset  (HP 

9872A) 

672 

IF  (PN.EQ.2)  THEN 

673 

PU-17 

1 

Plotter  Logical  Unit  =*  (HP 

7550A) 

674 

PXOFF-200 

1 

X-zero  Position  Offset  (HP 

7550A) 

675 

PYOFF-100 

1 

Y-zero  Position  Offset  (HP 

7550A) 

676 

END  IF 

677 

WIDTH-1,0 

f 

Character  Field  Width  (cm) 

678 

HEIGHT-2.0 

1 

Character  Field  Height  (cm) 

679 

Lmode-' SOLID’ 

! 

Solid  Line  Type  ON 

680 

Pattern-'  ' 

1 

Broken  Line  Type  OFF 

681 

Length-'  ' 

f 

Not  Used  For  Solid  Pattern 

682 

LF-CHAR(IO) 

f 

Linefeed  Character 

683 

CR-CHAR(13) 

! 

Carriage  Return  Character 

684 

C-'  ,  ' 

I 

Data  Separator  In  HP-CL  Commands 

685 

Etx-CHAR(3) 

1 

Labe!  Mode  Terminator  (HP  Default)' 

686 

Size-'N' 

1 

Plot  Size  Default  Is  NOTEBOOK 

687 

Vaxes-' 5' 

1 

Pen  Speed  (cm/sec)  For  Axes 

Vect  ors 

688 

Vspec— ' 5 ' 

! 

Pen  Speed  (cm/sec)  For  Spectrum 

689  C 
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690 

691 

692 

693 

694 
69^ 

696 

697 

698 

699 

700 

701 

702 

703 

704 

705 

706 

707 

708 

709 

710 

711 

712 

713 

714 

715 

716 

717 

718  - 

719 

720 

721 

722 

723 

724 

725 

726 

727 

728 

729 

730 

731 

732 

733 

734 

735 
■'36 

737 

738 

739 

740  C 

741 

742 


WRITE  (1,*)  DOWN , DOWN . DOWN 
470  WRITE  (1.480)  UP, ERASE Press  RETURN 
480  FORMAT  (T10,A2,A2,A31 ,A,A) 

READ  (1,99)  Icode 
!F  Hcode.NE.'  ')  CO  TO  47^ 

WRITE  (1. ' (T10,A2,A2.A21)’)  UP, ERASE 
WRITE  (PU,*, IOSTAT-N,ERR-9999)  ’DF;’ 
UXl-1800 
UX2-15500 
UYl-1200 
UY2-9500 

490  IF  (Size.EQ. ’F’ )  THEN 
Psize-'FULL' 

XP1-1350-PXOFF 
XP2-10000-PXOFF 
YPl-lOOO-PYOFF 
YP2-7500-PYOFF 
CO  TO  500 
END  IF 

Psize-' NOTEBOOK' 

XP1-1350-PXOFF 
XP2-8500-PXOFF 
YPl-lOOO-PYOFF 
YP2-7000-PYOFF 


if  plotter  is  ON:  ',BELL,’_' 


'Initializing  Plotter:' 


Default  User  Set  Plotting  Boundaries 
PI  -  (1800,1200)  &  P2  -  (15500,9500) 
For  B-size  Paper  On  HP  7550A  Plotter 


FULL  A-size  Paper  Plotting  Boundaries 
PI  -  (1350,1000)  &  P2  -  (10000,7500) 
Offset  P1,P2  For  HP  7550A  Plotter 


!  NOTEBOOK  Size  Plotting  Boundaries 
!  PI  -  (1350,1000)  &  P2  ~  (8500,7000) 
!  Offset  P1,P2  For  HP  7550A  Plotter 


LU  #33  is  device  #2  on  IEEE-488  interface  #1 

LU  #17  is  serial  device  on  HP  MUX  (9600  Baud,  XON/XOFF) 

500  WRITE  (PU,FMr-510) 

510  FORMAT  ( ' SPl ; PU; AP;TL1 ; DI 1 . 0 ; IW; ' ) 

WRITE  (PU.*)  ' IP' .XP1.C,YP1,C,XP2.C,YP2, ’ : ' 

Xleft-'O' 

Xright-' 10000' 

XOFF-XMIN 

IF  (XMIN.GT.XMAX)  THEN 
Xleft-' -10000' 

Xright-'O' 

XOFF-XMAX 
END  IF 

XMULT-T4/ ( XMAX-XM I N ) 

Ylow-'O' 

Ytop-' 10000' 

WRITE  (PU,*)  ' SC  .Xleft , 

RATIO-0.6 

LW I DTH-W I DTH*RAT 1 0 
LHE I CHT-HE I GHT*RAT I 0 
RAT 10-4. 0/3.0 
BW I DTH-W I DTH*RAT I 0 
BHE I CHT-HE I GHT*RAT 1 0 
WRITE  (PU,*)  'SR' .WIDTH, 


Scales  User  Units  To  0-10000  Range 
YMULT  Is  Set  In  Y-scale  Routines 
For  A  Scaled  Range  Of  0-10000  Units 
.Xright , ' , ' , Y1 ow , ' , ' , Yt op , ' ; ' 


LITTLE  Size  Characters  Default  To 
0.6  Of  The  User  Set  Character  Size 


!  BIG  Size  Characters  Are  Kept  In 
!  Proportion  To  The  User  Set  Size 
.HEIGHT, '  : VS'  , Vaxes .  '  :  ' 


520  WRITE  (1,*)  HQME,CLR,'_' 
CALL  Line(NCOL) 
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^43  WRITE  (1,530)  'CODE' , 'FUNCTION' , 'MIN' . 'MAX’ INC 

144  530  FORMAT  (T4 , A4 . T1 8 , A8 ,T40 . A3 ,T50 , A3 , T60 , A3) 

145  CALL  Line(NCOL) 

146  WRITE  (1,540)  'PL'.' . Plot  Limits . ',WN,WX 

M7  WRITE  (1,540)  'XS','..  ..X-scale . '  .  XM I N  ,  XMA.Y 

M8  IF  (Xcode.EQ. 'N' )  WRITE  ( 1 , ’ (T20 , A2 , A5) ' )  UP,'  (nm) ’ 

749  IF  (Xcode.EQ. 'W' )  WRITE  ( 1 , * (T20 , A2 , A14) ' )  UP,'  (cm-I)  x  lE-3' 

750  WRITE  (1,550)  'XT’,' . X-tick  Spacing . '  ,  FXT ,  LXT  ,  Xt  i  ck 

751  WRITE  (1,550)  'XL',’ . X-label  Spacing . ’  ,  FXL ,  LXL ,  XI  abe  ! 

752  WRITE  (1,*) 

753  WRITE  (1,560)  'YS',' . Y-scale . ', 

754  6c(-YOFF*YMULT)  .  (YSCALE-YOFF)*YMULT 

755  IF  (Ycode.EQ. 'E' )  WRITE  ( 1 , ’ (T20 , A2 . A9 , A5) ' )  UP,'  (Ext)  x  ',Mult 

756  WRITE  (1,570)  'YT'.' . Y-tick  Spacing . ', 

757  6t(FYT-YOFF)*YMULT,  (LYT-YOFF)*YMULT,  Yt  i ck*YMULT 

758  WRITE  (1.570)  'YL',* . Y-label  Spacing . ', 

759  &(FYL-YOFF)*YMULT, (LYL-YOFF)*YMULT,Ylabel*YMULT 

760  WRITE  (1.*) 

761  WRITE  (1.580)  'CS',' . Char. Size  (Width , He i ght )', WI DTH , HEIGHT 

762  WRITE  (1,600)  'PV',' . Pen  Velocity  (Axes , Spec) ', Vaxes , Vspec 

763  WRITE  (1,590)  'LT',' . Line  Type . '  ,  Lmode  ,  Pat  t  ern . 

764  £cLength 

765  WRITE  (1,590)  'AX',’ . Axes  Size . ',Psize,'  ','  ' 

766  WRITE  (1.610)  'AN',' . Annotate  Plot . ’ 

767  WRITE  (1,610)  'PA',' . PLOT  AXES . ' 

768  WRITE  (1,610)  'PS',’ . PLOT  SPECTRUM . ' 

769  WRITE  (1,610)  'DD',' . DIGITIZE/DISPLAY  DATA. ' 

770  WRITE  (1,610)  'EX’,' . EXIT  PLOTTING  MODE . ' 

771  540  FORMAT  (T5 , A2 ,T7 , A30 ,T40 , F4 . 1 ,T50 , F4 . 1 ) 

772  550  FORMAT  (T5 , A2 , T7 , A30 . T40 . F4 . 2 ,T50 . F4 . 2 , T60 . F4 . 2 ) 

773  560  FORMAT  (T5 , A2 , T7 , a30 , T40 . F5 . 3 , T50 . F5 . 3) 

774  570  FORMAT  (T5 , A2 , T7 , A30 , T40 , F5 . 3 , T50 , F5 . 3 , T60 , F4 , 3) 

775  580  FORMAT  (T5 , A2 , T7 . A30 .T40 , F4 . 2 ,T50 . F4 . 2) 

776  590  FORMAT  (T5 , A2 , T7 , A30 ,T40 , A8 , T50 , A6 , T60 . A6) 

777  600  FOF.’viAT  (T5  .  A2  ,  T7  ,  A30  ,T40 ,  A2  ,T50 ,  A2) 

778  610  FORMAT  (T5 , A2 ,T7 ,A30) 

779  CALL  Line(NCOL) 

780  WRITE  (1,*) 

781  620  WRITE  (1,*)  UP, ERASE, 

782  WRITE  (1 , ' (T3,A15.A,A2) ' )  'Enter  the  CODE BELL 

783  READ  (1,80)  Pcode 

784  CALL  Upper (Pcode) 

785  C  - 

786  IF  (Pcode. EQ. 'PL' )  CO  TO  700 

787  IF  (Pcode. EQ. 'XS' )  CO  TO  720 

788  IF  (Pcode. EQ. 'XT' )  GO  TO  730 

789  IF  (F-ode.EQ. 'XL' )  CO  TC  740 

790  IF  (Pcode. EQ. 'YS' )  CO  TO  760 

791  IF  (Pcode. EQ. 'YT' )  CO  TO  810 

792  IF  (Pcode. EQ. 'YL' )  CO  TO  820 

793  IF  (Pcode. EQ. 'CS' )  CO  TO  840 

794  IF  (Pcode. EQ. 'PV )  CO  TO  850 

795  IF  ( Pcode. EQ. 'LT' )  GO  TO  860 


796  IF  (Pcode.EQ. ' AX' )  GO  TO  890 

797  IF  (Pcode.EQ. 'EX’ )  CO  TO  980 

798  IF  (Pcode.EQ. 'PA' )  CO  TO  1000 

799  IF  (Pcode.EQ. 'AN' )  CO  TO  1100 

«0n  IF  (Pcode.EQ. 'PS')  CO  TO  2000 

801  IF  (Pcode.EQ. 'DD' )  CO  TO  3000 

802  CO  TO  500 

803  C  - 

804  700  WRITE  (1,*)  UP, ERASE,'  Plot  Limits:  (.Min, Max)  ',BELL.'_' 

805  READ  (1 ,*,ERR-700)  WN,WX 

806  CALL  Order(WN,WX) 

807  710  IF  ((WN.LT.WMIN) .OR. (WN.CT.WMAX))  WN-WMIN 

808  IF  ((WX.LT.WMIN) .OR. (WX.GT.WMAX))  WX-WMAX 

809  CO  TO  500 

810  C  - 

811  720  WRITE  (1,*)  UP, ERASE,’  X-scale;  (Min, Max)  ’,BELL,'_' 

812  READ  (1 ,*,ERR-720)  XMIN,XMAX 

813  XN— XMIN  !  X-axis  Is  Bidirectional  So  Only 

814  XX-XMAX  !  XN,XX  Are  Ordered  For  Tick  &  Label 

815  CALL  Order (XN, XX) 

816  CALL  Xaxis(XN,XX)  !  Set  Limits,  Tick  &  Label  Spacing 

817  GO  TO  710 

818  C  - 

819  730  WRITE  (1,*)  UP, ERASE,'  X-tick:  (Fi rst , Last , Space )  ',BELL,'_' 

820  READ  (1 ,*,ERR-730)  XN.XX.Xtick 

821  CALL  Order(XN,XX) 

822  IF  ((XN.LT.LMIN) .OR. (XX.GT.UMAX))  CO  TO  730 

823  Xtick-ABS(Xtick) 

824  •  FXT-XN 

825  LXT-XX 

826  Xlabel-Xtick 

827  CO  TO  750 

828  C  - - - 

829  740  WRITE  (1,*)  UP, ERASE, '  X-label:  (Fi rst , Last , Space)  ',BELL,'_' 

830  READ  (1 ,*,ERR-740)  XN.XX.XIabel 

831  CALL  Order(XN,XX) 

832  IF  ((XN.LT.LMIN) .OR. (XX. GT.LMAX))  CO  TO  740 

833  Xlabel-ABS(Xlabel ) 

834  750  FXL-XN 

835  LXL-XX 

836  GO  TO  500 

837  C  - 

838  760  WRITE  (1,*)  UP, ERASE,'  Y-scale:  (Min. .Max)  ',BELL.'_' 

839  READ  ( 1 , * , ERR-760)  YN.YX 

840  CALL  Order(YN,YX)  !  Y-scale  Is  Unidirectional 

84 1  EOFF— YN 

842  IF  (Ycode.NE. 'E' )  C'^  TO  800 

843  770  WRITE  (1,780)  UP. ERASE,'  Ext .Coeff. Scale ;  ',Mult, 

844  :  (U. . .UP  /  D. . .Down  /  A. . .Accept)  ?  ',BELL,'_’ 

S45  780  FOR’4AT  (A2  ,  A2  ,  A20  .  A6  ,  A37  ,  A .  A) 

846  READ  (1,99)  I  code 

847  CALL  Upper(Icode) 

848  IF  (Icode.EQ. 'A' )  GO  TO  800 


849 

850 

851 

852 
«5? 

854 

855 

856 

857 

858  790 

859 

860 
861 
862 

863  800 

864 

865 

866 

867  C 

868  810 

869 

870 

871 

872 

873 

874 

875 

876  C 

877  820 

878 

879 

880 

881  830 

882 

883 

884  C 

885  840 

886 

887 

888 

889  C 

890  850 

891 

892 

893 

894 

895 

896 

897 

898 

899  C 


IF  (Icode.EQ. 'U' )  THEN 
EMULT-10.0 
CO  TO  790 
END  IF 

IF  (  lrod<»  FO  '  O' )  THF'N 
EMULT-0 . 1 
CO  TO  790 
END  IF 
CO  TO  770 

ESCALE-Ei>CALE*EMULT 
RAT I  O-ESCALE/SCaLE 
YMULT-RAT I O/FACTOR ( J ) 

CALL  Exponent (RATI 0,Mul t ) 
CO  TO  770 

YSCALE- ( YX- YN ) /YMULT 
YOFF-EOFF/YMULT 
CALL  Yax is (YSCALE, Ycode) 
CO  TO  500 


!  Update  Ext .  Coeff .  Scale 
!  Power  Of  10  For  Exponent 
!  Scaling  Factor  For  Ext.  Coeff. 
!  Convert  Exponent  To  String 

!  YMULT  Alters  Scaling  For  The 
!  Extinction  Coefficient  Mod<» 

!  Set  Tick  &  Label  Spacing 


WRITE  (1,*)  UP, Erase,'  Y-tick;  (Fi rst , Last , Space) 
READ  (1 ,*,ERR-810)  YN,YX,Ytick 

CALL  Order(YN,YX) 

Yt  i  ck-ABS  ( Yt  i  ck)  AMULT 

FYT- ( YN+EOFF)  AMULT 

LYT-  ( YX+EOFF )  Al^ULT 

Ylabe 1-Yt ick 

CO  TO  830 

'  ,BELL. '_' 

WRITE  (1,*)  UP, ERASE,'  Y-IabeT;  (First.,  Last  ,  Space) 
READ  (1  ,*,ERR-820)  YN,YX,YUbel 

CALL  Order(YN,YX) 

Y 1  abe  1  -ABS  ( Y 1  abe  1 )  AT^ULT 

FYL- (YN+EOFF)  AI^ULT 

LYL- ( YX+EOFF )  A^ULT 

CO  TO  500 

' , BELL , ’ _ ’ 

WRITE  (1,*)  UP, ERASE,'  Char. size:  (Widt h , He i ght )  ' 
READ  (1 ,*,ERR-840)  WIDTH, HEICHT 

.BELL, '_' 

IF  (Pcode . EQ. ' AN' )  GO  TO  1140  !  Return  To  Ann^tat 

GO  TO  500 

ion  Mode 

WRITE  (1,*)  UP, ERASE,'  Pen  Velocity,  (Axes , Spect rum) :  ',BELL,'_ 
READ  (1 ,*,ERR-850)  X,Y 

IF  ((X.LT. 1 .0) .OR. (X.GT.36.0))  CO  TO  850 

IF  ((Y.LT. 1 .0) .OR. (Y.GT.36.0))  GO  TO  850 

CALL  Str(X, String, 2) 

Vaxes-St  r ing(2 : 3) 

CALL  Str(Y, String, 2) 

Vspec-S t  r i ng(2 ; 3 ) 

GO  TO  500 
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900 

901 

902 

903 

904 

905 

906 

907 

908 

909 

910 

911 

912 

913 

914 

915 

916 

917 

918 

919 

920 

921 

922 

923 

924 

925 

926 

927 

928 

929 

930 

931 

932 

933 

934 

935 

936 

937 

938 

939 

940 

941 

942 

943 

944 

945 

946 

947 

948 

949 

950 

951 

952 


860  WRITE  (1,*)  UP. ERASE. •  Line  Type:  (S . Solid,  B . Broken) 

&BELL. 

READ  (1.99)  I  code 
CALL  Uppe  r ( I code ) 

IF  ncode.EQ.  'S'  )  THEN 
Lmode-' SOLID' 

Pattern-'  ' 

Length-'  ' 

WRITE  (PU,*)  'LT;' 

GO  TO  500 
END  IF 

IF  (Icode.NE. 'B' )  GO  TO  860 
Lmode- 'PATTERN' 

870  WRITE  (1.*)  UP, ERASE.'  Pattern#:  (1-6)  ',BELL,'_' 

READ  (1 , ' (12) ' ,ERR-870)  N 
IF  ((N.LT. 1) .OR. (N.GT.6))  GO  TO  870 
NUMBER-FLOAT (N) 

CALL  Str(NUMBER, String. 1) 

Pattern(l : 1)-'#' 

Pattern(3:3)-String(2:2) 

880  WRITE  (1,*)  UP. ERASE,'  Pattern  Length:  (0-10)%  ',BELL,'_' 
READ  (1 ,*,ERR-880)  NUMBER 

IF  ((NUMBER. LT. 0.0) .OR. (NUMBER. GT. 10.0))  CO  TO  880 
CALL  Str(NUMBER, String, 2) 

K-2 

DO  WHILE  (String(K:K) .NE. '  ') 

K-K+1 
END  DO 

Le  ngt  h-S  t  r i ng ( 2 : K) 

WRITE  (PU,*)  'LT' ,Pattern(3:3) .Length, ': ' 

Length(K:K)-'%' 

GO  TO  500 


890  WRITE  (1,*)  UP, ERASE,'  F...FULL,  N ...  NOTEBOOK ,  U...USER  SET  ?  ’ 
■•DELL,  '_' 

READ  (1,99)  Size 
CALL  Upper (Size) 

IF  (Size.EQ. 'F' )  CO  TO  490 
IF  (Size.EQ. 'N' )  GO  TO  490 
IF  (Size.NE. 'U' )  CO  TO  890 
Psize-'USER  SET’ 

900  WRITE  (1,*)  UP, ERASE,'  Lower  Left:  (X,Y)  -  ' 

WRITE  (1,910)  UXl , ' , ' ,UY1 ,BELL, '_' 

910  FORMAT  (15, A, 15,'  :  A... Alter  or  RETURN  ?  '.A, A) 

READ  (1,99)  I  code 
CALL  Upper (I code) 
ir  (Icode.EQ. '  ' )  CO  TO  940 
IF  (Icode.NE. 'A' )  CO  TO  900 
920  WRITE  (1,930)  UP, ERASE, BELL. ' ?  ' 

930  FORMAT  (T23 , A2 , A2 , A , A3) 

READ  (1 ,*,ERR-920)  UXl.UYl 

UXl-ABS(UXl) 

im-ABS(UYl) 
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953  940  XPl-UXl 

954  YPl-UYl 

955  IF  (XPl.GT. 16000)  GO  TO  920 

956  IF  (YPl .GT. 11400)  GO  TO  920 

957  950  WKITE  O,*)  UP, ERASE,*  Upper  Right:  (X,Y)  -  ' 

958  WRITE  (1,910)  UX2, ' , ' ,UY2,BELL, 

959  READ  (1,99)  I  code 

960  CALL  Upper (I code) 

961  IF  (Icode.EQ.*  ')  GO  TO  970 

962  IF  (Icode.NE. 'A* )  GO  TO  950 

963  960  WRITE  (1,930)  UP, ERASE, BELL, ' ?  ' 

964  READ  (1 ,*,ERR-960)  UX2.UY2 

965  UX2-ABS(UX2) 

966  UY2-ABS(UY2) 

967  970  XP2-UX2 

968  YP2-UY2 

969  IF  (XP2.GT. 16000)  GO  TO  960 

970  IF  (YP2.GT. 10100)  GO  TO  960 

971  IF  ((XP1.GE.XP2).0R.(YP1.GE.YP2))  GO  TO  900 

972  IF  ((XP2-XP1) .LT. 1000)  GO  TO  900 

973  IF  ( (YP2 -YPl ) .LT. 1000)  GO  TO  900 

974  GO  TO  500 

975  C  - - - 

976  980  WRITE  (PU,*)  'SPO;' 

977  CLOSE  (PU, IOSTAT-N,ERR-9999) 

978  CO  TO  MENU 

979  C 

980  C  . . . - 

981  C 

982  C  AXES  PLOTTING  ROUTINE 

983  C 

984  C  -  Scaling  is  ON  - 

985  C 

986  C  Y-axis:  0-10000  user  units 

987  C  X-axis:  0-10000  user  units  (increasing  ->  right) 

988  C  -10000-0  user  units  (decreasing  ->  right) 

989  C 

990  C  -  Only  integral  position  values  are  sent  in  PA  commands 

991  C  for  compatibility  with  early  model  HP  plotters. 

992  C 

993  C  -  XMULT  &  YMULT  variables  scale  data  to  10000  digits 

994  C  for  full  scale,  ensuring  that  no  round-off  errors 

995  C  affect  plotting  resolution  on  early  model  HP  plotters 

996  C  which  do  not  accept  decimal  fractions  in  SC  or  PA. 

997  C 

998  C  -  Real  &  Integer  values  are  sent  within  the  HP-GL  commands 

999  C  as  appropriate  -  some  computers  send  extra  nulls  in  these 

1000  C  modes  causing  errors  on  early  HP  plotters.  In  such 

1001  C  cases  the  values  must  be  converted  to  string  literals 

1002  C  first.  The  present  form  works  correctly  on  the  HPIOOO. 

1003  C 

1004  C  - 

1005  C 
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1006  1000  WRITE  (1 ,*)  UP, ERASE. '  Plot t i ng  Axes :  ',BELL,'_' 

1007  WRITE  (PU,*)  'VS' .Vaxes, ’ ;LT:PA' .Xleft ,C,YIow, ’  ;PD;  ' 

1008  WRITE  (PU,*)  'PA' .Xright ,C.Ylow,C,Xright ,C,Ytop,C, Xleft ,C,Ytop,C, 

1009  &Xleft .C.Ylow, ' ;PU; ' 

1010  DX— 1E-5*ABS (XMAX)  !  Delta  X  &  Y  ensure  completion  of 

1011  DY“1 E-5*ABS (YSCALE)  !  DO  loops  with  fractional  steps 

1012  C  - - 

1013  DO  1010  V-FXT.LXT+DX,Xt ick 

1014  I-R0((V-XOFF)*XMULT) 

1015  WRITE  (PU,*)  'PA' , I , ’ , ' ,Ylow, • ;XT; ’ 

1016  1010  CONTINUE 

1017  C  - 

1018  DO  1020  V-FXL,LXL+DX.Xlabel 

1019  I-R0((V-XOFF)*XMULT) 

1020  WRITE  (PU,*)  'PA' , I , ’ , ' ,Ylow, ' ; * 

1021  X-R2(V) 

1022  CALL  Str(X, String, 4) 

1023  C 

1024  C  Find  #  of  digits  in  string,  ignoring  sign  Sc  trailing  blanks 

1025  C 

1026  K-2 

1027  DO  WHILE  (String(K:K) .NE. '  ') 

1028  K-K+1 

1029  END  DO 

1030  X-(FLOAT(K-2)-0.33)/2.0 

1031  WRITE  (PU,*)  'CP' ,-X, ' , -1 .LB' ,String(2;K-l ) ,Etx 

1032  1020  CONTINUE 

1033  C 

1034 

1035 

1036 

1037 

1038 

1039 

1040 

1041 

1042 

1043 

1044 

1045 

1046 

1047 

1048 

1049  C 


1050  1030 

1051 

1052 

1053  1040 

1054  C 

1055 

DO  1040  V-nrT,LYT+DY,Yt  ick 
I-R0(V*T4/YSCALE) 

WRITE  (PU,*)  'PA' .Xleft ' 
CONTINUE 

'  . 1 , '  ; YT ; ' 

DO  1050  V-FYL,LYL+DY,Y label 

1056 

Y-R3(V-YOFF) 

1057 

IF  (Ycode.EQ. 'E' )  THEN 

1058 

Y-R2 ( ( V-YOFF) *ESCALE/FACTOR ( J ) /SCALE ) 

1059 

END  IF 

IF  (Xcode.EQ. 'N' )  THEN 
Xst ring- 'Wavelength  (nm)' 

COFF-7.5-0.33 
END  IF 

IF  (Xcode.EQ. 'W' )  THEN 
Xstring—' Wavenumber  x  10' 

COFF-7.5+0.75-0.33 
END  IF 

String-'  5000' 

IF  (XMULT.lt. 0.0)  String-' -5000' 

WRITE  (PU.*)  'PA' ,String:i :5) , ' , ' ,Ylow, ' ;SR' .BWIDTH.C, BHEICHT, 
WRITE  (PU,*)  'CP' ,-COFF. ' ,-2;LB' .Xstring, Etx, 

&' SR' .WIDTH, ' , ' .HEIGHT, ' ; ' 

IF  (Xcode.EQ. 'N' )  GO  TO  1030 
WRITE  (PU,*)  'CPO, .5:LB-3' ,Etx 
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r 


1060  I-R0(V*T4/YSCALE) 

1061  WRITE  (PU,*)  'PA'  .Xleft , '  , ■  , I ,  ■  ;  ■ 

1062  CALL  Str(Y. String. 4) 

1063  IF  (String. EQ.'  0.0')  String-'  O' 

1064  C 

1065  C  Find  #  of  digits  in  string,  ignoring  s i en  &  trailing  blanks 

1066  C 

1067  K-2 

1068  DO  WHILE  (String(R:K) .NE. '  ') 

1069  K-K+1 

1070  END  DO 

1071  Y-FLOAT(K-l)+.33 

1072  WRITE  (PU,*)  'CP' ,-Y. ' .-.25;LB' ,String(l :K-1) .Etx 

1073  1050  CONTINUE 

1074  C  - - - 

1075  WRITE  (PU,*)  'PA' ,Xleft , ' ,5000;SR' .BWIDTH, ' , ' ,BHEICHT, 

1076  &' ;CP-5,0;DI0, 1 ; ' 

1077  IF  (Ycode.EQ. 'E' )  THEN 

1078  K-5 

1079  IF  (Mult(5:5) .EQ. '  ')  K-4 

1080  IF  (Mult(5:5) .EQ. ' . ' )  K-1 

1081  Ystring-'  x  '//MuIt(l:K) 

1082  COFF-(FLOAT(K)+13.5)/2.0-1 .5 

1083  IF  (K.EQ.l)  COFF-COFF-0 . 5 

1084  WRITE  (PU,*)  'CP' ,-COFF, ' ,0;CP-.33, .5;DR0,-1 ;LB3' ,Etx 

1085  WRITE  (PU,*)  'DR0,1:CP.33,-.5;LB' .Ystring(l:X+4), '  (M'.Etx 

1086  WRITE  (PU,*)  'CPO, .25;SR' .WIDTH, ',’ ,HE1CHT, ' ;LB-r  , Etx 

1087  WRITE  (PU,*)  'SR' .BWIDTH, ',' .BHEICHT, ’ ;CP.5,-.25;LBcm'  , Etx 

1088  WRITE.  (PU,*)’  'CPO,  .25:SR' , WIDTH, .HEIGHT,  ’  :LB-I ', Etx 

1089  WRITE  (PU,*)  ' SR' .BWIDTH, ',' .BHEICHT, ' 

1090  WRITE  (PU,*)  'CP0,-.25;LB)' ,Etx. ' ;DI1 ,0; ' 

1091  GO  TO  1060 

1092  END  IF 

1093  IF  (Ycode.EQ. 'A' )  K-10 

'094  IF  (Ycode.EQ. 'R' )  K-13 

1095  IF  (Ycode.EQ. 'T' )  K-14 

1096  COFF-FLOAT(K)/2.0-0. 167 

1097  WRITE  (PU,*)  'CP' ,-COFF, ' ,0;LB' ,Ystring(l :K) , Etx, ' ;DI1 ,0; ' 

1098  1060  WRITE  (1 ,*) 

1099  IF  (YOFF.EQ.0.0)  GO  TO  1090 

1100  C  - 

1101  1070  WRITE  (1,*)  UP, ERASE,'  Plot  a  dashed  baseline,  (Y  or  N)  ? 

1102  &BELL,'_' 

1103  READ  (1,99)  I  code 

1104  CALL  Upper (I code) 

1105  IF  (Icode.EQ. 'N' )  GO  TO  1090 

1106  IF  (Icode.NE. 'Y')  CO  TO  1070 

1107  1080  WRITE  (1,*)  UP, ERASE,'  Dash  Length:  (1-5V',  •,BELL,’_' 

1108  READ  (1 ,*)  X 

1109  IF  ((X.LT. 1.0) .OR. (X.CT.5.0))  GO  TO  1080 

1110  IX-RO(X) 

nil  IZ-RO(YOFF*T4/YSCALE) 
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1112  WRITE  (PU,*)  'LT2.  MX,  '  ;PA'  .Xright  ,  '  ,  MZ.  '  ;PD;PA'  .Xleft  , 

1113  , IZ, ' ;PU;LT; ' 

1114  1090  IF  (Psize.EQ. 'NOTEBOOK')  CO  TO  1200 

1115  CO  TO  2080 

1116  C 

1117  C  - 

1118  C 

1119  C  Annotacion:  -  Labels,  Parameters,  Title 

1120  C 

1121  C  - 

1122  C 

1123  1100  WRITE  (1 ,*)  UP, ERASE,'  L.. .Labels,  P ...  Parameters ,  T.. .Title 

1124  ,  X. . .Exit  ?  ' .BELL, 

1125  READ  (1,99)  Icode 

1126  CALL  Upper (Icode) 

1127  IF  (Icode. EQ. 'T' )  CO  TO  1110 

1128  IF  (Icode. EQ. 'L' )  GO  TO  1120 

1129  IF  ( Icode. EQ. 'P' )  CO  TO  1200 

1130  IF  ( Icode. EQ. 'X' )  GO  TO  2080 

1131  GO  TO  1100 

1132  C  - - - 

1133  1110  TITLE-LABEL(J) 

1134  WRITE  (1 , ' (A2,A2,A4,A72,T76,A) ■ )  UP, ERASE,'  ?  <', TITLE, '>’ 

1135  WRITE  (1 , ' (T5,A2,A,A) ' )  UP,BELL,'_' 

1136  READ  (1,120)  TITLE 

1137  IF  (TITLE. EQ.'  ')  TITLE-LABEL( J ) 

1138  WRITE  (1,*)  UP, ERASE,'  Plot t ing  Ti t 1 e BELL 

1139  K-72 

1140  DO  WHILE  (TITLE(K:K) .EQ. '  ’) 

1141  K-K-1 

1142  END  DO 

1143  COFF-FLOAT(K)/2.0 

1144  IX-5000 

1145  IF  (XMULT.lt.  0.0)  IX— 5000 

1146  'V-1 

1147  IF  (Psize.EQ. 'NOTEBOOK' )  IY=2 

1148  WRITE  (PU,*)  ' PA'  , IX,C,Ytop, '  ;SR1 , 2;CP'  , -COFF.  '  ,  '  .  lY,  '  :  ' 

1149  WRITE  (PU,*)  'LB' ,TITLE(1 :K) ,Etx 

1150  CO  TO  1100 

1151  C  - 

1152  1120  WRITE  ( 1 , ' ( A2 , A2 , A1 3 ,T58 , A) ' )  UP, ERASE,'  Label:  ? 

1153  WRITE  (! , ' (T14,a2,A,A) ’ )  UP.BELL,'_’ 

1154  READ  (I , ' (A40) ' )  Text 

1155  IF  (Text.EQ.'  ')  CO  TO  1100 

1156  K-40 

1157  DO  WHILE  (Text (K : K) . EQ . '  ') 

1158  K  K-1 

1159  END  DO 

1160  1130  WRITE  (1  ,*)  UP, ERASE,  '  Label  ;  '  .Text  1  :K:  .  '  :  OK  ^Y  or  N  ■?  '  . 

1161  &BELL,'_' 

1162  READ  (1,99)  Icode 

1163  CALL  Upper (Icode) 

1164  IF  ( Icode. EQ. 'N' )  CO  TO  1120 

1165  IF  ( Icode .NE. 'Y' )  GO  TO  1130 


1166 

1167 

1168 

1169 

1170 

1171 

1172 

1173 

1174 

1175 

1176 

1177 

1178 

1179 

1180 
1181 
1182 

1183 

1184 

1185 

1186 

1187 

1188 

1189 

1190 

1191 

1192 

1193 

1194 

1195 

1196 

1197 

1198 

1199 

1200 
1201 
1202 

1203 

1204 

1205 

1206 

1207 

1208 

1209 

1210 
1211 
1212 

1213 

1214 

1215 

1216 

1217 

1218 


1140  WRITE  (1,*)  UP, ERASE,'  C. . . 
c'V. . .Vplot ,  X. . .Exit  ?  •  ,BE 
READ  (1,99)  Icode 
CALL  Upper (Icode) 

IF  ( Icode. EQ. 'C )  CO  TO  1150 
IF  ( Icode. EQ. 'M' )  GO  TO  1160 
IF  (Icode. EQ. 'V )  GO  TO  1170 
IF  (Icode. EQ. 'P' )  GO  TO  1180 
IF  (Icode. NE. 'X' )  CO  TO  1140 
GO  TO  1100 


C. . .Char .size 
' ,BELL, 


. Move , 


1150  WRITE  (1,*)  UP, ERASE,'  Char. Size  (Wi dt h , He i ght )  ?  ’,BELL,'_' 
READ  (1 ,*,ERR-1150)  LWIDTH, LHEICHT 
GO  TO  1140 

1160  WRITE  (1,*)  UP, ERASE,'  Coordinate,  0-100%:  (X,Y)  ?  '.BELL,'_' 
READ  (1 ,*,ERR-1160)  X,Y 

I X- RO ( X* 1 0 0 . 0+ ( XM I N -XOFF) *XMULT ) 

IY-R0(Y*100.0) 

WRITE  (PU,*)  'PA' , IX, ' , • , lY, ' ; • 

GO  TO  1140 

Special  Entry  Point  From  Digitize  Routine  (Pcode . EQ . ' DD' ) 
Special  Offset  Controls  Are  Used  For  Digitize  Mode  Labels 

1170  WRITE  (PU,*)  'DI0,1;'  !  Rotate  Labelling  Axis  90  Degrees 

1180  WRITE  (1,*)  UP, ERASE,'  C. . .Centered,  L...Left  Justified, 

&'R... Right  Justified  ?  ',BELL,'_' 

READ  (1,99)  Icode  • 

CALL  Upper (Icode) 

IF  ( Icode.  EQ.  'C  )  THEN 
X— FLOAT  (K)/2.0 
IF  (Pcode. EQ. 'DD' )  THEN 

X-X-0.5  !  Alter  Centering  For  Leading  Blank 

IF  (K.GT.8)  THEN  !  Alter  Centering  For  Labelling  Both 

1-2  !  X  &  Y  Coordinates  In  Digitize  Mode 

DO  WHILE  (Text ( I : I) .NE. '  ') 

I-I  +  l 

END  DO  !  Find  Space  Between  X  &  Y  Labels 

X— FLOAT  (I +  1) 

END  IF 
END  IF 
GO  TO  1190 
END  IF 

IF  ( Icode. EQ. 'R' )  THEN 
X— FLOAT  (K) 

IF  (Pcode.  EQ.  ’  D."' )  THEN 

X-X-2.0  !  Alter  Right  Justification  To  Stop 

END  IF  !  Peak  Label  Writing  Over  Spectrum 

GO  TO  1190 
END  IF 

IF  (Icode. NE. 'L' )  CO  TO  1180 

X- 0.0  !  Convert  X=*0 . 0  To  '0.0'  With  Str 


!  Find  Space  Between  X  &  Y  Labels 
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1 


1219 

1220 
1221 

1222  C 

1223  1190 

1224 

1225 

1226 

1227 

1228 

1229 

1230 

1231 

1232 

1233 

1234  C 

1235  1200 

1236  1210 

1237 

1238 

1239 

1240 

1241 

1242 

1243 

1244 

1245 

1246 

1247 

1248 

1249 

1250 

1251 

1252 

1253 

1254 

1255 

1256 

1257 

1258 

1259 

1260 

1261  C 

1262  C 

1263  C 

1264  C 

1265  C 

1266  C 

1267  C 

1268  C 

1269  C 

1270  C 

1271  C 

1272  C 


IF  (Pcode.EQ. 'DD')  THEN 

X“1.0  !  Alter  Left  Justification  To  Stop 

END  IF  !  Peak  Label  Writing  Over  Spectrum 


CALL  Str(X, String, 4) 

WRITE  (1,*)  UP, ERASE,’  Plotting  Label: '.BELL 

WRITE  (PU,*)  'SR' .LWIDTH, ' , ■ .LHEIGHT, ' ;CP' .Stringd :5) , ' .-.25; 

WRITE  (PU,*)  'LB' ,Text(l;k),Etx. ' ;CP0, .25:DI1,0; ' 

IF  (Pcode.EQ. 'DD' )  THEN 

I X-N I NT (T4*FLOAT ( I X-Xl ) /FLOAT (X2 -XI ) + ( XM I N-XOFF ) *XMULT ) 
IY-NlNT(T4*FLOAT( lY-Yl ) /FLOAT (Y2-Y1 ) ) 

WRITE  (PU.*)  'PA' , IX. ’ . ' , lY. ' ;PU; ' 

CO  TO  3000  !  Return  To  Digitize  Routine 

END  IF 
CO  TO  1100 


IF  (Psize.EQ. 'FULL')  CO  TO  1100 

WRITE  (1,*)  UP, ERASE,’  Plot  Parameters,  (Y  or  N)  ?  ',BELL,'_' 
READ  (1,99)  I code 
CALL  Upper (I code) 

IF  (Icode.EQ. 'N' )  GO  TO  1100 
IF  (Icode.NE. 'Y' )  CO  TO  1210 

WRITE  (1,*)  UP, ERASE,'  Plotting  Parameters BELL 
WRITE  (PU,*)  'PA' .Xright .C.Ytop, ' : ' 

WRITE  (PU,*)  'SR.75,2;CP0,-.5;LB' ,LF,LF, '  File:  ',Fname(J), 
&CR,LF,LF,Etx 

WRITE  (PU,*)  'LB  Date:  ' . DATE( J) , CR, LF, LF, Et x 
WRITE  (PU,*)  'LB  Cone. (M) : ' .CONC(J) ,CR,LF,LF,Etx 
WRITE  (PU,*)  'LB  Path  (cm) :  ’  .  PATH(  J)  ,  CR,  LF ,  LF‘,  Et  x 
WRITE  (PU,*)  'LB  Rate  (nm/sec):  ’, 

&Pstr(3,PARAM(J ,3)+l) ,CR,LF,LF.Etx 
WRITE  (PU,*)  'LB  Period  (sec):  ', 

&Pstr(15,PARAM(J, 15)+1) ,CR,LF,LF,Etx 
WRITE  (PU,*)  'LB  Ref.  Mode:  ', 

&Pstr(5,PARAM(J,5)+l) .CR.LF.LF.r  x 
WRITE  (PU,*)  'LB  Beam  Mode:  ', 

&Pstr(16,PARAM(J , 16)+I ) .CR, LF, LF, Etx 
WRITE  (PU,*)  'LB  SBW  (nm)  ; ' , VARI ABLE( J , 1 0) , CR , LF, LF , Et x 
WRITE  (PU,*)  'LB  Slit  Cain: ’ , VARI ABLE( J , 6) , CR , LF , LF , Et x 
WRITE  (PU,*)  'LB  Slit  Height;  ', 

£.Pstr(23,PARAM(J,23)+l)  .Etx 
CO  TO  1100 

END  OF  AXES  PLOTTING 
SPECTRUM  PLOTTING  ROUTINE 

-  Velocity  Select  -  5  cm/sec  fPaper) ,  10  cm/sec  (Acetate) 

-  Real  variables  are  used  in  'VS'  commands 

-  Scaling  in  ON  with  Absorbance  x  1E4  to  match  axes 

-  Both  Integer  and  String  Numbers  are  used  in  'PA'  commands 

-  HP  I/O  Subsystem  Performs  Binary  <=—>  ASCII  Conversions 
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1273  C 

1274  2000  M-1 

1275  WRITE  (1,*)  UP, ERASE,'  Increase  step  size,  (Y  or  N)  ?  ',BELL,'_' 

1276  READ  (1,99)  Icode 

1277  CALL  Upper (Icode)  < 

1278  IF  (Icode. EQ. 'N' )  CO  TO  2020 

1279  IF  (Icode. NE. 'Y' )  GO  TO  2000 

1280  2010  WRITE  (1,*)  UP, ERASE,'  Step  Mul t ip  I ier ,  (>-l):  ',BELL,'_' 

1281  READ  (1 , ' (13) ' ,ERR-2010)  M 

1282  IF  (M.LT.l)  M-1 

1283  2020  X-STEP*FLOAT(M)  i 

1284  CALL  Str(X, String, 4) 

1285  C 

1286  C  Find  #  of  digits,  ignoring  sign  &  trailing  blanks 

1287  C 

1288  K~2  j 

1289  DO  WHILE  (String(K;K) .NE. '  ')  ^ 

1290  K-K+1 

1291  END  DO 

1292  WRITE  (1,*)  UP, ERASE,'  Plot t ing  Spect rum  (3  ' , St r i ng (2 ; X-1 ) , 

1293  &'  nm  steps;  ',BELL,'_' 

1294  C  I 

1295  C  - - - - - 

1296  C 

1297  C  WINDOW  SETTING  ROUTINE 

1298  C 

1299  C  -  Soft  clip  limits  are  set  to  PI  &  P2  scaling  coordinates 

1300  C  so  that  off  scale  data  do  not  plot  on  graph.  This  should  I 

1301  C  not  happen  normally  since  tne  plot  boundaries  are  made  to 

1302  C  match  the  X-scale  limits  automatically  if  the  latter  are 

1303  C  smaller  than  the  scan  range. 

1304  C  -  Plotter  sends  output  parameters  as  A, B,C, D,CR, LF 

1305  C  where  (A-D)  are  string  integers  (ASCII),  which  are  read 

1306  C  into  integer  variables  (ASCII  translation  is  automatic).  I 

1307  C  -  Some  computers  handshake  on  CR  leaving  Li  in  buffer. 

1308  C  In  these  cases  use  READ  A,B,C,D,E  (N.B.  CHARACTER  E*l) 

1309  C  to  clear  plotter  buffer  before  the  next  READ  cycle. 

1310  C 

1311  C  - 

1312  C  ► 

1313  2030  WRITE  (PU,*)  ' VS ' ,Vspec , ' ;OP; ' 

1314  READ  (PU,*)  X1,Y1,X2,Y2 

1315  WRITE  (PU,*)  ' IW' ,X1 ,C,Y1 ,C,X2,C,Y2, ' ; ’ 

1316  IF  (Xcode.EQ. 'N' )  THEN 

1317  START-R0((WMAX-WX)/STEP)+1 

1318  FINISH-R0(FLOAT(ND)-(WN-WMIN)/STEP)  K 

1319  GO  TO  2040 

1320  END  IF 

1321  START-R0((T4/WMIN-T4/WN)/STEP)+I 

1322  FINISH-RO(FLOAT(ND)-(T4/WX-T4/WMAX)/STEP)  i 

1323  C  -  J 

1324  2040  I-O  '-1 

1325  RATI0-FACT0R(1)/FACT0R(2)  , 
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I 


I 


[ 


I 


1326 

1327 

1328 

1329 

C 

1330 

1331 

1332 

1333 

C 

1334 

1335 

2050 

1336 

1337 

1338 

1339 

1340 

2060 

1341 

1342 

2070 

1343 

C 

1344 

2080 

1345 

1346 

1347 

1348 

1349 

C 

1350 

C 

1351 

C 

1352 

C 

1353 

C 

1354 

C 

1355 

C 

1356 

C 

1357 

C 

1358 

C 

1359 

C 

1360 

c 

1361 

3000 

1362 

1363 

1364 

1365 

1366 

1367 

1368 

1369 

1370 

1371 

1372 

1373 

1374 

1375 

3010 

1376 

3020 

1377 

3030 

1378 

DO  2070  K-START, FINISH, M 

X-W(J,K)  !  Acquired  Wavelength 

IF  (Xcode.EQ. 'N' )  THEN 

X— WX-FLOAT( I )*STEP  !  Calculated  Wavelength 

1 X-RO ( ( X-XOFF ) *XMULT ) 

GO  TO  2050 
END  IF 

X-T4/WN-FLOAT(I)*STEP  !  Calculated  Wavelength 

I X-RO ( (T4/X-XOFF) *XMULT) 

IF  (Code.EQ. ’DS')  THEN 

IY-R0(((A(1 .K)-A(2.K)*RATIO)+YOFF)*T4/YSCALE) 

CO  TO  2060 
END  IF 

IY-R0((A(J  ,k:)+YOFF)*T4/YSCALE) 

WRITE  (PU,*)  ’ PA' , IX,C, lY. ■ :PD: ' 

I-I+M 

CONTINUE 


WRITE  (PU.*)  'PU;LT; ' 
Lmode-' SOLID' 
Pattern-'  ' 

Length-'  ' 

GO  TO  500 


PLOTTER  DIGITIZING  ROUTINE 

-  HP-CL  command  'OA'  provides  pen  position  without  using 
the  clumsy  ENTER  key  on  the  plotter. 

-  The  'DP'  &  'OD'  commands  are  used  with  ENTER  in  cases 
where  the  plotter  is  remote  from  the  keyboard. 


WRITE  (1,*)  HOME,CLR,'_' 

TITLE— ' Digit ize  Pen  Position' 

CALL  Center (TITLE) 

CALL  Line(NCOL) 

CALL  Di gi t ize (X, Xcode , Y, Ext , K) 

WRITE  (1,3010)  'X-position-  ' ,X,Xord 
String—'  Absorbance  -  ' 

IF  (Ycode.EQ. 'R' )  String-'  %  Reflect.  -  ' 
IF  (Ycode.EQ. 'T' )  String-'  %  Transmit.-  ’ 
WRITE  (1,3020)  String.Y 

IF  ((Ycode.EQ. 'A'). OR. (Ycode.EQ. ’E'))  THEN 
WRITE  (1,3030)  'Ext.Coeff.  -  ' 

WRITE  (1,*)  Ext(l;K).’  /Wcm' 

END  IF 

FORMAT  (/,T5,A14,T20,F5.3,A7) 

FORMAT  (/,T5,A14,T20,F6.4) 

FORMAT  (/,T5,A14) 

WRITE  (1,*)  DOWN 
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1379 

1380 

1381 

1382 

1383 

1384 

1385 

1386 

1387 

1388 

1389 

1390 

1391 

1392 

1393 

1394 

1395 

1396 

1397 

1398 

1399 

1400 

1401 

1402 

1403 

1404 

1405 

1406 

1407 

1408 

1409 

1410 

1411 

1412 

1413 

1414 

1415 

1416 

1417 

1418 

1419 

1420 

1421 

1422 

1423 

1424 

1425 

1426 

1427 

1428 

1429 

1430 

1431 

1432 


WRITE  (1,3040)  'D . Digitize  Position' 

WRITE  (1,3040)  'L . Label  At  Position' 

WRITE  (1,3040)  'S . Segment  Display  ' 

WRITE  (1,3040)  'X . EXIT  To  Plotting  ' 

3040  FORMAT  (/,T16,A23) 

WRITE  (1,*)  DOWN, DOWN 

3050  WRITE  (1 , ’ (T10,A2,A2,A6,A,A) • )  UP, ERASE, ' Code ; 
READ  (1,99)  I code 
CALL  Upper (I code) 

IF  (Icode.EQ. 'X’)  CO  TO  500 
IF  (Icode.EQ. 'D' )  CO  TO  3000 
IF  (Icode.EQ. ’S')  CO  TO  4000 
IF  (Icode.NE. 'L' )  CO  TO  3050 
3060  WRITE  (1,*)  UP, ERASE,'  Place  PEN  In  Stall  s^l  , 
&’Then  Press  RETURN  ',BELL,*_’ 

READ  (1,99)  I code 
IF  (Icode.NE.’  ')  CO  TO  3060 
WRITE  (PU,*)  'PU;SP1;' 

CALL  Digi t i2e(X,Xcode ,Y,Ext ,K) 

3070  WRITE  (1,*)  UP, ERASE,’  X...X-label,  Y. . .Y-lab 
&'  Q. . .Quit  ?  ' ,BELL, 

READ  (1,99)  I code 
CALL  Upper (I code) 

IF  (Icode.EQ. 'Q' )  CO  TO  3000 
IF  (Icode.EQ. 'Y* )  CO  TO  3080 

IF  ((Icode.NE. 'X' ) .AND. (Icode.NE. 'B' ))  CO  TO  3070 
CALL  Str(X, String, 5) 

Text-String  !  String*14,  Text*40 

K-40‘  !  Are  Reuseable  Strin 

DO  WHILE  (Text (K:K) .EQ. '  ') 

K-K-1  !  Find  Start  Of  Trail 

END  DO 

IF  (Icode.EQ. 'X' )  CO  TO  3110 


,BELL, 


Y.  . .Y-label 


. Both ' 


!  String*14,  Text*40  And  T1TLE*72 
!  Are  Reuseable  String  Variables 

!  Find  Start  Of  Trailing  Blanks 


3080  IF  (Ycode.EQ. 'E' )  THEN 
St  r i ng-Ext 
CO  TO  3090 
END  IF 

CALL  Str(Y, String, 5) 

3090  IF  (Icode.EQ. 'B' )  THEN 

TITLE-Text(l ;K)//'  (’//String 
Text-T1TLE(1  ■■40) 

GO  TO  3100 
END  IF 
Text-St  r i ng 
3100  K-40 

DO  WHILE  (Text(K;K) .EQ. ’  ’) 

K-K-1  !  Find  Start  Of  Trailing  Blanks 

END  DO 

IF  (Icode.EQ. 'B' )  THEN 
Text(K+l :K+2)-'  )' 

K-K+2 
END  IF 

3110  GO  TO  1170  !  Transfer  To  Vertical  Label  Routine 
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1433 

1434 

1435 

1436 

1437 

1438 

1439 

1440 

1441 

1442 

1443 

1444 

1445 

1446 

1447 

1448 

1449 

1450 

1451 

1452 

1453 

1454 

1455 

1456 

1457 

1458 

1459 

1460 

1461 

1462 

1463 

1464 

1465 

1466 

1467 

1468 

1469 

1470 

1471 

1472 

1473 

1474 

1475 

1476 

1477 

1478 

1479 

1480 

1481 

1482 

1483 

1484 

1485 

1486 


C 

c 

c 

c 

c 


Display  Segment 


Of  Spectrum 


4000  String-'  (nm) ' 

IF  (Xcode.EQ. ’W)  String-'  (cm-1)  x  lE-3’ 

WRITE  (1,*)  UP, ERASE,'  Scan  Range;  ',WM1N,'  to',WMAX,'_' 

WRITE  (1,*)  String, DOWN, DOWN 

4010  WRITE  (1,*)  UP, ERASE,'  DISPLAY  Range:  (MIN, MAX)  ?  ',BELL,'_' 
READ  (1 ,*,ERR-4010)  W1,W2 
CALL  Order(Wl,W2) 

IF  ((Wl .LT.WMIN) .OR. (W1 .GT.WMAX))  CO  TO  4010 
IF  ( (W2. LT.WMIN) .OR. (W2.CT.»flMAX))  GO  TO  4010 
IF  (Xcode.EQ. 'N' )  THEN 

START-R0( (WMAX-W2)/STEP)+1 
FINISH-RO(FLOAT(ND)-(W1-WMIN)/STEP) 

CO  TO  4020 
END  IF 

START-RO ( ( T4/WM IN  T  4/Wl ) /STEP)  + 1 

FINIS  H-RO ( FLOAT ( ND ) - ( T4/W2 -T4 /WMAX) / STEP ) 

4020  IF  ((Ycode.EQ. 'R' ) .OR. vVcode.EQ. 'T' ))  THEN 
Text-Ystring 

IF  (Ycode.EQ. 'R' )  Text-'  '//Ystring 

TITLE-' Energy  Wavelength  ' //Text ( 1 : 15) 

CO  TO  4030 
END  IF 

TITLE-' Energy  Wavelength  Absorbance  Ext.Coeff.  (/M/cm)' 
4030  WRITE  (1,*)  HOME,CLR.'_' 

CALL  Center (TITLE) 

CALL  Line (60) 

L-1 

DO  4050  I-START, FINISH 
X-T4/W(J, I) 

Y-A(J, I) 

IF  (Code.EQ. 'DS' )  Y-A( 1 . I ) -A(2 , 1 ) *FACTOR( 1 ) /FACTOR(2) 

WRITE  (1,*) 

IF  ((Ycode.EQ. 'A' ) .OR. (Ycode.EQ. 'E' ) )  THEN 
WRITE  (1,4060)  X,W(J, I) ,Y,Y/FACTOR(J) 

ELSE 

WRITE  (1,4070)  X,W(J. I) ,Y 
END  IF 
L-L+1 

IF  (L.EQ. 10)  THEN 
WRITE  (1 ,*)  DOWN 

4040  WRITE  (1,4080)  UP, ERASE, ' Press  RETURN  to  Continue  ',BELL.'_' 
READ  (1,99)  I  code 
IF  (Icode.NE.'  ')  GO  TO  4040 
WRITE  (1,*)  HOME,CLR,'_' 

CALL  Center(TITLE) 

CALL  Line(60) 

L-1 
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1487  END  IF 

1488  4050  CONTINUE 

1489  4060  FORMAT  (T9 , F7 . 3 ,T21 , F7 . 2 .T35 , F7 . 4 , T49 , CIO . 4) 

1490  4070  FORMAT  (T14 , F7 . 3 , T30 . F7 . 2 . T47 , F7 . 2) 

1491  4080  FORMAT  (TIO , A2 , A2 , A25 . A, A) 

1492  WRITE  (1.*)  DOWN 

1493  4090  WRITE  (1,4080)  UP. ERASE. ' Press  RETURN  for  MENU  ',BELL,'_' 

1494  READ  (1.99)  I  code 

1495  IF  (Icode.NE.'  ')  CO  TO  4090 

1496  CO  TO  3000 

1497  C 

1498  C  - 

1499  C 

1500  C  Exit  Program 

1501  C 

1502  C  - 

1503  C 

1504  9000  WRITE  (1.*)  UP, ERASE, UP 

1505  STOP 

1506  C 

1507  C  - 

1508  C 

1509  C  IEEE-488  Error  Exit 

1510  C 

1511  C  — . . . . . 

1512  C 

1513  9999  WRITE  (1,*)  '  Error  #' .N 

1514  STOP 

1515  END 

1516  C 

1517  C  ***********************  end  of  MAIN  PROGRAM  ************************* 

1518  C 

1519  C  Print  a  TITLE  Centered  in  72  columns 

1520  C 

i  521  C  ************************************************************************** 

1522  C 

1523  SUBROUTINE  Center (TITLE) 

1524  INTEGER  I,J,N 

1525  CHARACTER  TITLE*72 , BLANK*36 

1526  BLANK-' 

1527  1-72 

1528  J-0 

1529  DO  WHILE  ( ICHAR(TITLE( I ; I ) ) . EQ. 32) 

1530  J-J+1 

1531  I-72-J 

1532  END  DO 

1533  N-J/2 

1534  WRITE  (1,*)  BLANK(1 :N) ,T1TLE(1 : 1) 

1535  RETURN 

1536  END 
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1537  C 

1538  C 

1539  C 

1540  C 

1541  C 

1542  C 

1543  C 

1544 

1545 

1546 

1547 

1548 

1549 
15:50 

1551 

1552 

1553 

1554 

1555 

1556 

1557 

1558 

1559 

1560 

1561 

1562 

1563 

1564 

1565 

1566 

1567 

1568 

1569 

1570 

1571 

1572 

1573 

1574 

1575 

1576 

1577 

1578 

1579 

1580 

1581 

1582 

1583 

1584 

1585 

1586 
587 


Extinction  Coefficient  Rounding 


SUBROUTINE  Coeff (E, String, K) 

INTEGER  F.K.M 
REAL  E,S 

CHARACTER* (*) S t r i ng 
CHARACTER  Mu It *3 
F-0 
S-1 .0 

IF  (E.LT.0.0)  S— 1.0 
E-ABS(E) 

IF  ((E.NE.0.0) .AND. (E.NE.1.0))  GOTO  10 
String-'  0.0000' 

IF  (E.EQ.1.0)  String-'  1.0000' 

K-7 

GO  TO  50 

10  IF  ((E.GT. 1 .0) .AND. (E.LT. 10.0))  THEN 
K-K+l 
GO  TO  30 
END  IF 

TF  (E.GE.10.0)  GO  TO  20 
DO  WHILE  (E.LT. 1 .0) 

E-E*10.0 
F-F-1 
END  DO 
GO  TO  30 

20  DO  WHILE  (E.CE. 10.0) 

E-E/10.0 
F-F+1 
END  DO 
30  E-E*S 

CALL  Str(E,String,K) 

40  K-2 

DO  WHILE  (String(K:K) .NE. '  ') 

K-K+l 
END  DO 

IF  (F.EQ.O)  GO  TO  50 
E-FLOAT (F) 

CALL  Str(E,Mult ,2) 

IF  (Mult (1 : 1) .EQ. '  ’)  Mul t ( 1 : 1 )-' +' 

M-3 

IF  (Mult(3:3) .EQ. '  ' )  M-2 
String(K;K+M)-'E'//Mult(l :M) 

K-K+M 
50  RETURN 
END 


Digitize  Pen  Position 


1588  C 

1589  C 

1590  C 

1591  C 

1592  C 

1593  C 

1594  C 

1595 

1596 

1597 

1598 

1599 

1600 
1601 
1602 

1603 

1604 

1605 

1606 

1607 

1608 

1609 

1610 
1611 
1612 

1613 

1614 

1615 

1616  C 

1617  C 

1618  C 

1619  C 

1620  C 

1621  C 

1622  C 

1623 

1624 

1625 

1626 

1627 

1628 

1629 

1630 

1631 

1632 

1633 

1634 

1635 

1636 

1637 

1638 

1639 

1640 


SUBROUTINE  Di gi t ize (X.Xcode , Y, Ext .K) 

INTEGER  IX. lY, 12, J.K.PU.Xl ,X2,Y1 ,Y2 

REAL  FACTOR(2) , NUMBER, XMIN.XMAX.YOFF.YSCALE.X.Y 

CHARACTER*(*)  Ext ,Xcode 


CO\WON  /DIGIT/FACTOR,XMIN,XMAX,YOFF,YSCALE, IX, lY, J ,PU,X1 ,X2,Yl ,Y2 
WRITE  (PU,*)  'OP; ' 

READ  (PU,*)  X1.Y1.X2.Y2 
WRITE  (PU,*)  'OA;  ' 

READ  (PU,*)  IX.IY.IZ 

X- ( XMAX-XM I N ) *FLOAT ( I X-Xl ) /FLOAT ( X2 -XI ) +XM I N 
IF  (Xcode.EQ. 'W' )  THEN 
X-ANINT(X*1000.0) 

GO  TO  10 
END  IF 

X-ANINT(X*10.0)/10.0 

10  Y-YSCALE*FLOAT( lY-Yl ) /FLOAT(Y2-Yl ) -YOFF 

K-4  !  K  -  #  Of  Digits  Precision 

NUMBER-Y/FACTOR ( J ) 

CALL  Coe ff (NUMBER, Ext , K)  !  K  -  #  Of  Characters  In  String 

RETURN 

END 


Extinction  Coefficient  Scale  Multiplier 


SUBROUTINE  Exponent (N ,Mul t ) 
INTEGER  I 
REAL  M,N 

CHARACTER  String*14 
CHARACTER* (*)  Mult 
M-0.0 
N-ABS(N) 

I-INT(N+.5) 

IF  (I.EQ.l)  THEN 
Mult-'l 
GO  TO  20 
FND  IF 

IF  ( 1 . GT . 1 )  CO  TO  10 
DO  WHILE  ( I .LT. 1) 

N-N*10.0 
M-M-1 .0 
I-INT(N+.5) 

END  DO 


I 

I 


■ 


I 


•J 


•J 


59 
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1641  CALL  Str(M, String, 2) 

1642  Mult-'lE-'//String(2:3) 

1643  CO  TO  20 

1644  10  DO  WHILE  (I .CT. 1) 

1645  N-N/10.0 

1646  M-M+1.0 

1647  I-INT(N+.5) 

1648  END  DO 

1649  CALL  Str(M, String, 2) 

1650  Mult-' IE+’//String(2;3) 

1651  20  RETURN 

1652  END 

1653  C 

1654  C  - 

1655  C 

1656  C  Print  a  line  of  N  characters  (72  columns  max) 

1657  C 

1658  C  - 

1659  C 

1660  SUBROUTINE  Line(N) 

1661  INTEGER  I.N 

1662  CHARACTER  BLANK*72 , DLINE*72 . SPACE*36 

1553  SFA5E— ' 

1664  BLANK-SPACE//SPACE 

1665  SPACE-' - ' 

1666  DLINE-SPACE//SPACE 

1667  IF  (N.GT.72)  N-72 

1668  I-(72-N)/2 

1669  WRITE  (1,*)  BLANK(1 : n .DLINE(1 :N) 

1670  RETURN 

1671  END 

1672  C  ’ 

1673  C  - - — - - 

1674  C 

1675  C  Order  MIN, .MX  Entries  If  Required 

1676  C 

1677  C  - 

1678  C 

1679  SUBROUTINE  Order (MIN .MAX) 

1680  REAL  MIN, MAX, SWAP 

1681  IF  (MIN. LT. MAX)  CO  TO  10 

1682  SWAP-MIN 

1683  MIN-MAX 

1684  MAX-SWAP 

1685  10  RETURN 


Rescale  Data  Array  Between  Absorbance  And  Transmission 


1687  C 

1688  C 

1689  C 

1690  C 

1691  C 

1692  C 

1693  C 

1694  $EMA/DATA/ 

1695  SUBROUTINE  Rescale(J ,K,ND) 

1696  INTEGER  J.K.ND 

1697  REAL  A(2, 10001) 

1698  CONWON  /DATA/A 

1699  IF  (K.EQ.l)  GO  TO  20 

1700  DO  10  I-l .ND 

1701  A(J, I)-ALOG10(100.0/A(J. 1)) 

1702  10  CONTINUE 

1703  RETURN 

1704  20  DO  30  I-l ,ND 

1705  A(J, I)-100.0/10.0**(A(J, I)) 

1706  30  CONTINUE 

1707  RETURN 

1708  END 

1709  C 

1710  c  - : - 

1711  C 

1712  C  Convert  Number  To  ASCII  String 

1713  C 

1714  C  . . - . . . . 

1715  C 

1716  SUBROUTINE  St r (VALUE . S t r i ng, PREC) 

1717  INTEGER  ASCI  I , DECPT , 1 . J , LENSTR.NDI GIT , NUMBER , PREC 

1718  REAL  VALUE 

1719  DOUBLE  PRECISION  DEC I MAL , FRACT I  ON , TEN 

1720  CHARACTER  Concat*14 , Digi t ( 12) , S i gn , St r i ng*l 4 

1721  LOGICAL  INTEGER 

1722  C  LOGICAL  TEST 

1723  INTEGER-. TRUE. 

1724  C  TEST-. FALSE. 

1725  DECPT-0 

1726  J-0 

1727  TEN-10.0 

1728  Sign-'  ' 

1729  Concat-'  ' 

1730  C  IF  (TEST)  WRITE  (1,*)  ’  Value  Entered  =  ’.VALUE 

1731  IF  (VALUE.lt. 0.0)  Sign-'-' 

1732  IF  (VALUE. EQ. 0.0)  GO  TO  100 

1733  DEC I MAL- ABS (VALUE) 

1734  DO  WHILE  (DECIMAL . GE . 1 . 0) 

1735  DECIMAL-DECIMAL/TEN 

1736  J-J+1 

1737  END  DO 

1738  DECPT- J 

1739  C  IF  (TESTl  WRITE  (1.*)  ’  »  of  Whole  Digits:  '.DECPT 


1740  IF  (DECPT.EQ.O)  GO  TO  30 

1741  DO  20  J-1,DECPT 

1742  DECIMAL-DECIMAL*TEN 

1743  NUMBER- I NT (DECIMAL) 

1744  ASCn-NUMBER+48 

1745  Digit(J)-CHAR(ASCI I) 

1746  FRACTION-DECIMAL-NUMBER 

1747  DEC I MAL-D I NT ( FRACT I ON*TEN** ( PREC- J ) + . 5 ) /TEN** ( PREC - J ) 

1748  20  CONTINUE 

1749  C  IF  (.NOT. TEST)  GO  TO  30 

1750  C  WRITE  (1,*)  '  The  Whole  Digits  -  ',(Digit(I),  I=I,DECPT) 

1751  30  J-DECPT 

1752  C  IF  (TEST)  WRITE  (1,*)  '  Decimal  Fraction  -  '.DECIMAL 

1753  IF  (DECIMAL. NE. 0.0)  INTEGER-. FALSE. 

1754  IF  (DECPT.GE. 12)  GO  TO  40 

1755  DO  WHILE  (DECIMAL. NL. 0 . 0) 

1756  J-J+1 

1757  DECIMAL-DECIMAL*TEN 

1758  NUMBER- I NT (DECIMAL) 

1759  ASCI I-NUMBER+48 

1760  Digit(J)-CHAR(ASCII) 

1761  FRACTION-DECIMAL-NUMBER 

1762  DECIMAL-D I NT ( FRACT I ON*TEN** ( PREC- J ) + . 5 ) /TEN** ( PREC- J ) 

1763  IF  (DECIMAL. EQ. 1 .0)  THEN 

1764  DIGIT(J)-CHAR(ASCI I+l) 

1765  DECIMAL-0.0 

1766  END  IF 

1767  IF  (J.GE.12)  DECIMAL-0.0 

1768  END  DO 

1769  40  NDICIT-J 

1770  C  IF  (.NOT. TEST)  CO  TO  50 

1771  C  WRITE  (1,*)  '  The  Characters  -  ’,(Digit(I),  I-1,NDICIT) 

177^  50  IF  (NDICIT.CT. 12)  CO  TO  200 

1773  DO  60  I-l .NDICIT 

1774  Concatd  :  I)-Digit(I) 

1775  60  CONTINUE 

1776  IF  (INTEGER)  GO  TO  80 

1777  IF  (DECPT.EQ.O)  CO  TO  70 

1778  String-Sign//Concat (1 iDECPT)//' . ' //Concat ( DEcPT+l : 14) 

1779  RETURN 

1780  70  String-S ign//' . ' //Concat 

1781  RETURN 

1782  80  String-Sign//Concat 

1783  RETURN 

1784  100  String-'  0.0' 

1785  RETURN 

1786  20G  WRITE  (1,*)  '  Error  in  data:  (too  many  digits)' 

1787  STOP 

1788  END 


Convert  String  Entry  To  Uppercase  If  Required 


1789  C 

1790  C 

1791  C 

1792  C 

1793  C 

1794  C 

1795  C 

1796  SUBROUTINE  Upper(Code) 

1797  INTEGER  LENSTR.N 

1798  CHARACTER*(*)  Code 

1799  LENSTR-LEN(Code) 

1800  DO  10  I-l .LENSTR 

1801  N-lCHAR(Code(I  •  D) 

1802  IF  (N.CT.96)  Code( I ; I )-CHAR(N-32) 

1803  10  CONTINUE 

1804  RETURN 


Convert  ASCII  String  To  Numeric  Value  (10  Digits  Max'm) 


1806  C 

1807  C 

1808  C 

1809  C 

1810  C 

1811  C 

1812  C 

1813 

1814 

1815 

1816 

1817 

1818 

1819 

1820  C 

1821 

1822  C 

1823 

1824 

1825 

1826 

1827 

1828 

1829 

1830  C 

1831  C 

1832 

1833 

1834 

1835 

1836 

1837 

1838 

1839 

1840 

1841 

1842 

1843 

1844 

1845 

1846 

1847 

1848 

1849 

1850 

1851  C 

1852 

1853 


SUBROUTINE  Val ( St r i ng , VALUE) 

INTEGER  DECPT,EXPON,LENSTR.N,NUM(10) 

REAL  VALUE 

DOUBLE  PRECISION  MULT ,  S I GN .  TEN ,  DEC  1 NIAL 

CHARACTER  Ascii 

CHARACTER*(*)  String 

LOGICAL  INTEGER 

LOGICAL  TEST 

INTEGER-. TRUE. 

TEST-. FALSE. 

J-1 

K-0 

DECPT-0 

SIGN-1.0 

TEN-IO.O 

DECIMAL-0.0 

LENSTR-LEN(String) 

IF  (TEST)  IVRITE  (1,*)  '  String  Number  -  '.String 
IF  (TEST)  WRITE  (1,*)  '  String  Length  -'.LENSTR 
DO  100  I-l, LENSTR 
Asci i-St  r ing( I ; I ) 

N-lCHAR(Asci i) 

IF  ((N.GE.48) .AND. (N.LE.57))  CO  TO  20 
IF  (N.EQ.46)  INTEGER-. FALSE. 

IF  (N.EQ.46)  DECPT-K 
IF  (N.EQ.45)  SIGN— 1.0 
CO  TO  100 
20  NUM(J)-N-48 
K-J 
J-J  +  1 

100  CONTINUE 

IF  ((DECPT.EQ.O) .AND. (INTEGER),  DECPT-K 
DO  200  J-1 ,K 
EXPON-DECPT-J 
MULT-TEN**EXPON 
DEC  I MAL-DEC I  NiAL+NUM  ( J  )  '■VMULT 
200  CONTINUE 

VALUE-S  I  GN*DEC  I  .MAL 

IF  (TEST)  WRITE  (1,*)  '  Value  -'.VALUE 

RETURN 

END 
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1854  C 

1855  C 

1856  C 

1857  C 

1858  C 

1859  C 

1860  C 

1861 
1862 

1863 

1864 

1865 

1866 

1867 

1868 

1869 

1870  C 

1871 

1872 

1873 

1874 

1875 

1876 

1877 

1878 

1879 

1880  C 

1881 
1882 

1883 

1884 

1885 

1886  C 

1887 

1888 

1889 

1890 

1891 

1892 

1893 

1894 

1895 

1896 

1897 

1898 


Set  X-axis  Plotting  Parameters 


SUBROUTINE  Xaxis(XL,XH) 

REAL  LMIN,LMAX,WN,WX,Xt ick , FXT, LXT, XI abe 1 , FXL , LXL , XH , XL , XS ( I  7 ) 
COMMON  /XPARAM/LMIN.LMAX,\W,WX,Xt ick,FXT,LXT,Xlabel , FXL, LXL 
DATA  (XS(I) . I-l , 17)/.05. . 1 , .2, .25. .5, 1 . ,2. ,2.5,5. , 10. .20. ,25. , 
&50. , 100. ,200. ,250. ,500./ 

LMIN—XL  !  Set  Limits  To  X-scale  Range 

LMAX— XH  !  Defined  By  XL,XH  Arguments 

WN-XL 
WX-XH 


Xt ick-0. 1*(XH-XL)  !  Set  Xtick  To  1/10  X-scale 

IF  ((Xtick.LT.XS(2)) .OR. (Xtick. CT.XS(17)))  THEN 

FXT-XL+Xtick  !  Test  For  Xtick  Outside  The 

LXT-XH-Xtick  !  Preferred  Interval  Range 

Xlabel— Xt ick 

FXI,-FXT 

LXL-LXT 

CO  TO  10  !  Use  Default  1/10  Interval 

END  IF 


I-l 

DO  WHILE  (Xt ick. GT . XS ( I ) )  !  Compare  1/10  Scale  Xtick  Value 

I-I+l  !  To  Find  The  Nearest  Preferred 

END  DO  !  Interval  <  1/10  Of  X-scale 

Xt ick-XS( I-l) 


FXT-Xt ick*AINT(XL/Xt ick)+Xt ick  ! 
LXT-Xt ick*AINT(XH/Xt ick)  ! 

IF  (ABS(XH-LXT) .LT.0.01)  THEN 
LXT-LXT-Xt ick  ! 

END  IF 

XI abe 1-Xt i ck*2 . 0  ! 

FXL-X 1 abe I *A 1 NT ( XL/X 1 abe 1 ) +X 1 abe I 
LXL-Xlabel*AINT(XH/Xlabe 1 ) 

IF  (ABS(XH-LXL) .LT.0.01 )  THEN 
LXL-LXL-X label 
END  IF 
10  RETURN 


Truncate  ->  FXT  <=■  XL 

Truncate  ->  LXT  <=■  XH 

Decrement  I f  LXT  -  XH 

Truncate  ->  FXL  <=  XL 

Truncate  ->  LXL  <=■  XH 


r+Xt ick) 


('  +X1  abe  I 


!  Decrement  I f  LXL  =  XH 


Set  Y-axis  Plotting  Parameters 


SUBROUTINE  Yaxi s (YSCALE . Ycode) 

REAL  Ytick.FYT.LYT.Ylabel .FYL.LYL, YSCALE, YS(17) 

CHARACTER  Ycode 

COMMON  /YPARAMAt  ick.FYT.LYT,  Y  label  ,  FYL.LYL 

DATA  (YS(I) . I-l , 17)/. 0005, . 001 , . 002 , . 005 , . 01 , . 02 , . 05 , . 1 , . 2 , . 5 , 

&1.  ,2.  ,5.  ,10. ,20. ,50. ,100./ 

Yt ick-0. 1*YSCALE 

IF  ((Yt  ick.LT.YS(2))  .OR.  (Yt  ick.  CT .  YS  ( 17) )  .OR.  (Ycode. EQ.  'E'  ))  THE.N 
FYT-Yt i ck 
LYT-YSCALE-Yt ick 
Ylabel-Yt ick 
FYL-FYT 
LYL-LYT 
GO  TO  10 
END  IF 


DO  WHILE  (Ytick.CT.YS(l))  ! 

I-I4l  ! 

END  DO  ! 

Yt ick-YS(I-l) 

FYT-Yt ick 

LYT-Yt ick*AINT(YSCALE/Yt ick) 

IF  (ABS(YSCALE-LYT) .LT. lE-4)  THEN 
LYT-LYT-Yt ick 
END  IF 

Ylabel-Yt ick*2 . 0 
FYL-Y label 

LYL-Y 1 abe 1 *A I NT ( YSCALE/Y label) 

IF  (ABS(YSCALE-LYL) .LT. lE-4)  THEN 
LYL-LYL-Y label 
END  IF 
RETURN 
END 


Compare  1/10  Scale  Ytick  Value 
To  Find  The  Nearest  Preferred 
Interval  <  1/10  YSCALE 


Truncate  ->  LYT  <=  YSCALE 
Decrement  If  LYT  -  YS 


Truncate  ->  LYL  <=  YSCALE 
Decrement  If  LYL  =  YSCALE 
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Wait  Specified  Delay  (sec) 


1941  C 

1942  C 

1943  C 

1944  C 

1945  C 

1946  C 

1947  C 

1948  SUBROUTINE  Wait (DELAY) 

1949  REAL  DELAY, PERIOD .Tzero .Time 

1950  PERIOD-0.0 

1951  Tzero-Time( I )  | 

1952  DO  WHILE  (PERIOD. LT. DELAY) 

1953  PERIOD-Time(I)-Tzero 

1954  END  DO 

1955  RETURN 

1956  END 

1957  C 

1958  C  - 

1959  C 

1960  C  Read  Time  (sec)  from  the  HP  1000' s  RTE-6  Operating  Svstem 

1961  C 

1962  C  Note:  I  is  a  dummy  argument,  no  values  are  passed 

1963  C 

1964  C  - 

1965  C 

1966  REAL  FUNCTION  Time(I) 

1967  INTEGER  ICODE, ITIME(5) 

1968  I CODE-11 

1969  CALL  EXEC( ICODE, ITIME) 

1970  Time-FLOAT(ITIME(1))/I00.0+FLOAT(1T1ME(2))+FLOAT(ITIME(3))-'^60.0 

1971  &+FLOAT(lTIME(4))*3600.0 

1972  RETURN 

1973  END 
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Appendix 
Sample  Plots 


i 
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LD  'i: 


(■^_LU3  f-N)  —  ^  J 
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Example  Of  Resca 


